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Abstract 

We discuss the physics case of a high luminosity i3-Factory running at the T(5S') 
resonance. We show that the coherence of the B meson pairs is preserved at this 
resonance, and that B s can be well distinguished from and charged B mesons. 
These facts allow to cover the physics program of a traditional i?-Factory and, at 
the same time, to perform complementary measurements which are not accessible 
at the T(45). In particular we show how, despite the experimental limitations 
in performing time-dependent measurements of B s decays, the same experimental 
information can be extracted, in several cases, from the determination of time- 
integrated observables. In addition, a few examples of the potentiality in measuring 
rare B s decays are given. Finally, we discuss how the study of B s meson will improve 
the constraints on New Physics parameters in the B s sector, in the context of the 
generalized Unitarity Triangle analysis. 
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1 Introduction 



The study of B mesons at the 5-Factories has considerably improved our understanding of 
flavour physics in the Standard Model (SM), confirming that CP violation is well described 
by the Cabibbo-Kobayashi-Maskawa (CKM) matrix [1]. In particular, the agreement 
between the measured value of sin 2/3 in b — > ccs decays [2, 4] and the prediction from 
the indirect constraints on the CKM parameters p and fj [3] can be considered the first 
precision test of the CKM mechanism. Additional constraints have been added, leading to 
a remarkable improvement of the determination of the parameters of the CKM matrix [5]. 

A large set of these measurements is based on the study of the coherent B-B state, 
produced by the decay of the T(45). The possibility of studying simultaneously time- 
dependent CP asymmetries and rates in several B decays to charged and neutral particles 
pushed the physics program of BaBar [6] and Belle [7] well beyond the initial expectations. 
The key ingredients of the success of BaBar and Belle are the large achieved luminosities 
and the clean environment in which B mesons are reconstructed. The two collaborations 
are expected to improve the situation by doubling their datasets in the next two years. 

Recently, the study of B physics received an additional boost by the results from 
the Tevatron experiments CDF [8] and D0 [9]. In particular, the measurement of the 
oscillation frequency of the B s -B s system, Am s [10] and the comparison to the prediction 
from the indirect determination [4] represents an additional test of the SM, having a value 
comparable to the measurement of sin 2/3. In the next future, new measurements will be 
added, thanks to the increased dataset of the Tevatron experiments and the start of the 
LHCb program [11]. 

Since these measurements will be performed at hadronic machines, where Bd, B u , B s 
and B c mesons are simultaneously produced, it will be possible to cover a much larger 
physics sample. In particular, from the study of B s mesons it is possible to extract some of 
the fundamental quantities that are also accessible at the 5-Factories (as the CKM phase 
7 or New Physics (NP) parameters) with reduced theoretical uncertainty with respect to 
the case of Bd mesons. In addition, thanks to the different quark content of the initial 
state, several B s decays, which are comparable to interesting decay modes of the Bd 
meson, can provide more experimental information. For instance, the decay B s — > KK 
has the same interesting features of Bd — > Kir, but since it is a CP eigenstate, a study of 
the time-dependent CP asymmetry can be performed. Another example is B s — ► J /ipcf), 
similar to Bd — > J/ipK*°, but which can be reconstructed in a CP eigenstate with all 
charged particles in the final state and therefore with a higher efficiency. 

However, we stress that in these latter cases the measurements will be performed in a 
very different context with respect to the clean e + e~ environment. Because of this, several 
measurements that can be accessed at the 5-Factories, such as those involving neutral 
particles (i.e. 7r°'s, rf decays, radiative photons, etc.) can not be carried out with the same 
high accuracy as in an e + e~ facility. We thus suggest here to investigate the possibility 
of performing them at an e + e~ environment. Such opportunity might be provided by 
a next generation B- Factory, running also at the T(5S) resonance [13]. We will show 

1 The study of — > J / tpK*° in CP eigenstate requires the K*° mesons to be reconstructed from ^5^° 
combinations, with a loss of about a factor of six in efficiency and a large background, because of the 7r° 
in the final state. Due to the low reconstruction efficiency even at the B-Factories the time-dependent 
study of this channel [12] is affected by a large statistical error. 
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Figure 1: CLEO measurement of R = o~(e + e — > qq)/o~(e + e — > /i/x) as a function of the 
invariant mass W, in the region ofW ~ 11 GeV [16]. 



how a run at the T(5S) could provide, with respect to already approved experiments, 
both complementary informations and independent determinations of interesting CKM 
parameters and quantities sensitive to NP. 

In the recent past, CLEO and Belle performed short runs at the T(5S), measuring 
the main features of this resonance with a few collected data. These first results give just 
a flavour of the potentiality of a B- Factory running at the T(5S), which is the subject of 
this paper. 

We assume the same detector performances of the existing 5-Factory experiments and 
we use a full detector simulation to study the physics potential of running at T(55), as a 
function of the integrated luminosity. 

We show that it is possible to separate B s mesons from and charged B mesons, 
even in presence of several photons in the final state. In this way, running at the T(5S), 
it would be possible to perform the same physics measurements than at the present B- 
Factories, after imposing kinematic cuts to reject B s events. 

2 The T(5S) Production and Decays 

The T(5S) resonance is a J CP = 1 state of a bb quark pair, having an invariant 
mass mr(5s) = (10.865 ± 0.0008) GeV [14]. The cross section of T(5S) production in 
e + e~ collisions is a(e + e~ — > T(5S)) = 0.301 ± 0.002 ± 0.039 nb [15]. This resonance was 
discovered by CLEO at the CESR e + e~ collider, measuring the total cross section above 
the T(4S) resonance [16]. The distribution of R = <r(e + e~ — > qq)/a(e + e~ — > fifi) as a 
function of the invariant mass is shown in Fig. [TJ 

The knowledge on the properties of this bound state comes from the 0.42 fb _1 collected 
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Table 1: BR ofT(5S) decays, measured by CLEO [15] and BELLE [17] collaborations. 
The last column shows the values used throughout this paper. 



by the CLEO experiment at CESR [16] and by the 1.86 fb _1 collected by the Belle 
detector at KEKB, during the T(55) engineering run [17]. The Belle collaboration has 
also collected a sample of about 20 fb™ 1 during June 2006. 

Unlike the T(4S") state, this resonance is heavy enough to decay in several B meson 
states. In particular, it can decay to vector- vector (VV), pseudoscalar- vector (PV), and 
pseudoscalar-pseudoscalar (PP) combinations of i?**^, B d *^ and Bi*^ mesons. The T(5<S') 
resonance can also decay into B^B^n. 

In the case of PP decays, the two B mesons are produced in a 1 state, as in the 
case of the T(4S") resonance. When the produced mesons are neutral, they exhibit flavour 
oscillations, because of the coherence of the initial state, the oscillation frequency being 
determined by the mass differences Am dyS . 

The case of PV and VV decay is more complicated. The vector states B* (q = u, d, s) 
decay in B* B^ through an electromagnetic interaction which, up to a very good level 
of approximation, can be considered as an instantaneous process. Under this assumption, 
it can be shown (see Sec. 14. ip that the time evolution of VV states is similar to that of 
the PP state. For events generated by PV decays, a difference comes from the opposite 
eigenvalue of the charge-conjugation operator. 

Quark models [18] predict the production to be mainly in the excited final states 
B ds B ds with B ds — >• B d;S j. Even if relative branching ratios (BR) of T(5S) are not 
precisely known yet, the current measurements confirm this picture. We summarize in 
Tab.[T]the corresponding experimental results, as reported by CLEO [15] and BELLE [17], 
along with the values used in this paper. For unmeasured values, we made educated 
guesses based on the predicted (and observed) VV dominance. 

The numbers quoted in the table, together with the different values of the cross 
sections, show that B d and charged B mesons can be produced at the T(5S'), with a 
rate about six times smaller than at the T(4S). This fact should be kept in mind: a 
reduction in the statistics of B d and charged B mesons is the price to pay in order to 
study B s particles at a £>-Factory. 
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Figure 2: Distribution of L 2 /L for B s — > J/if>(f> events simulated at the T(5S) resonance 
in a VV (dots) and PP (black line) configurations. 



3 Reconstruction of B—B pairs 

The reconstruction of B-B pairs at the T(5S) proceeds like at the traditional B- 
Factories. The only relevant difference is the fact that several final states can be produced, 
with different momenta. 

In the case of B s mesons, the largest fraction (~ 94% of the events) is produced in 
a VV mode. One can then tune the event reconstruction and selection on these events, 
considering the PV and the PP modes as a source of background. To avoid decreasing 
the reconstruction efficiency, one should not look for either the photons produced by the 
B* d decay or the one or more pions produced in the continuum bb production. 

At a traditional 5-Factory, BB events are distinguished from qq background (q = 
u, d, s, c) using a set of variables related to the distribution of the decay products in the 
center-of-mass (CM) system of the T resonance. We consider here the two quantities L 
and L 2 , defined as [19] 

L,- = J>*||cos0*P (7 = 0,2), (1) 

i 

where p* is the momentum of particle i in the e + e~ rest frame, 9* is the angle between 
p* and the thrust axis of the B candidate and the sum runs over all reconstructed par- 
ticles except for the S-candidate daughters. In data analyses, Lq and L 2 are typically 
combined into a Fisher discriminant T [20]. Alternatively, the ratio L 2 /Lq is used [21]. 
The discriminating power against the qq background for these two variables is very close, 
so that we use one or the other in the subsequent sections. 

At the T(5S), Lq and L 2 preserve all their features, even in presence of additional 
unreconstructed particles (photons and pions) produced with the BB pair. This is shown 
in Fig. El where the distribution of the ratio L 2 /Lq for B s — > J/ip(j) generated events is 
given for PP (continuous line) and VV (dots) events. We do not distinguish between VV, 
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PV, and PP in the following, when dealing with topological variables. 

Kinematic variables are also used to reject background. At the T(4S), a BB event is 
usually characterized by the two following variables: 

• AE = E* B — y/s/2, the energy difference between the reconstructed B candidate in 
the center of mass (CM) frame E* B and one half the total CM energy y/s. 

• the beam-energy substituted mass titles = \/ (s/2 + pi • p-B) 2 /Ef — pb 2 , where the 
B momentum pb and the four-momentum (Ei,pi) of the e + e~ initial state are 
defined in the laboratory frame. 

These two variables are found to be Gaussian distributed and almost uncorrelated for 
those B decays having only charged tracks in the final state. For B mesons coming from 
T(4S) decays, ities (AE) peaks at the value of the B d mass (at zero). The presence 
of photons in the final state (coming from decays of intermediate particles, radiative B 
decay or bremsstrahlung) introduces a correlation (typically less than 20%), induced by 
the missing energy. In these cases, it was found [21] that a better signal identification is 
achieved replacing ities with m miss , defined as m miss = \p e + e - — Pb\, where p e + e - is the 
four- momentum of the e + e~ initial state and Pb is the four-momentum of the B candidate 
after applying a mass constraint on it. 

To illustrate the kinematic properties of B mesons coming from T(5S) decays, we 
consider the case of T(55) — ► B*B* events. After the two B* — ► B s ^ decays, the final 
state is represented by the decay products of the two B s mesons and by the two additional 
photons which are not reconstructed. When calculating AE from the observed decay 
products, E* B assumes the value of the energy of the B s meson, while the value to use for 
a/s/2 corresponds to half the energy of the B*B* system. This difference generates a shift 
of the AE distribution of —47 MeV, corresponding to the mass difference -mg.. In 
addition, the B s momentum distribution is smeared by the spread of the photon energy 
in the CM frame, since the photon is monochromatic only in the B* rest frame, while 
the laboratory frame is boosted @ and the energy of the photon becomes a function of 
the polar angle in the CM frame. At the same time, m,Es values are shifted by the same 
quantity to larger values. This is shown in Fig. [31 where the distribution of Bd, s — > J/ip^ 
events in the AE vs tties plane is given. For B s decays the three bumps (corresponding 
to PP, PV and VV decays) are aligned along a straight line, identified by the relation: 
AE = —(rriEs — m B s )- Similar considerations are valid for Bj and charged B mesons. 

For m m i SS the situation is more complicated, since the calculation of this variable 
requires a mass hypothesis for the reconstructed B (B^ or B s meson). One is forced to 
define two variables, w^ iss and rn s miss , and to use one or the other according to the case 
under study. 

In the following, we present three examples of B decays at the T(5S) resonance. We 
start from the simple case of a final state without neutral particles (Bd. s —>■ J/^<p). After 
showing the good performances in separating VV events from PV and PP events and B s 
from Bd, we consider the case of B final states with two photons (B^ s — ► iv Kg) and 
three photons (B^ s —>■ K* 0, y with K*° —>■ 7t°Kg). In principle, one would expect these 
decays (and in general all the decays to final states with photons) to be more problematic, 

2 In order to be able to study also Bd physics, T(55) decays should be studied at an asymmetric 
B-Factory, where also time-dependent measurements of Bd decays can be performed. 
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Figure 3: Distribution of AE vs. rriEs for a sample of simulated B ds — > J/ipcf) decays at 
the T(55) resonance. Events coming from Bq*^ B q *^ (q — d,s) are all generated with the 
same relative rate. We use full boxes for q = d and empty boxes for q = s. The color scale 
identifies VV, VP and PP events (from the darker to the lighter). Events from B d B d -n 
events are also shown (full dark boxes on the bottom-right) . 



since only a part of the energy of the photons is detected in the calorimeter (introducing 
asymmetries in the shapes of the kinematic variables). As we show, the use of m m i SS 
instead of tues can provide a good separation even in this case. Other choices adopted 
in literature (such the use of rather than AE) do not give advantages and are not 
considered in the following. 

3.1 Selection of Bd, s — ► J/if><fi 

Fig. [3] shows the distribution of signal Monte Carlo events for B ds — > J/ip<fi decays, 
obtained assuming all the bb decays of the T(5S) to have the same production rate. 
The values of the two variables are calculated from the daughters of the reconstructed 
B mesons, without searching for additional particles (photons or pions) produced in the 
T(55) decay tree. 

It can be noted that the PP, PV and VV configurations are well separated in the 
(rriEs,AE) plane. B ds B ds events are centered at AE ~ while B ds B ds and B^ s B* ds are 
shifted to negative values of AE and higher values of tties due to the missing energy of 
the photons originating from B* — > Bgj decays. In a similar way, BBtt events are further 
shifted. Three important features should be noted: 

• The resolution of both variables is worse in the case of B d mesons than for B s 
mesons, since B d mesons receive a larger momentum than B s mesons in the CM 
frame of the T(5S), resulting in a wider distribution of the B momentum in the 
laboratory frame. 

• The resolution for VV events is worse than for PP ones, since the presence of the 
photon introduces an additional source of energy spread. 
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Figure 4: Top row : distribution of ities for a sample of simulated B<i (left) and B s 
(right) decaying to J/ip<p at the T(5S) resonance. Bottom row : distribution of m d miss for 
a sample of simulated B^ (left) and m s miss for a sample of B s (right) decaying to J /ip<p at 
the T(5S) resonance. From the left to right, we show events coming from B q B q , B q B*, 
B*B q (dashed) and B*B* (q = d,s), all generated with the same relative rate. In writing 

BqBq the first B q *^ is intended as the reconstructed one. 




Figure 5: Distribution of tues fl e ft) an d AE (right) for Bd — > J/ip<p candidates, from 
a sample of simulated T(55) — > BdBdF events (line). B s — > J/ip<fr events, generated in 
T(5S) — * B*B* events (dots), are also shown for comparison. 



• In the case of PV events, the distribution is the sum of two components: a broad 
one, due to reconstructed B mesons originating from the B* decay; and a narrow 
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Figure 6: Top row : distribution of m^s for a sample of simulated B d (left) and B s 
(right) decaying to K^tt at the T(5S) resonance. Bottom row : distribution of m d miss for 
a sample of simulated B& (left) and m s miss for a sample of B s (right) decaying to K^tt at 
the T(5S) resonance. From the left to right, we show events coming from B q B q , B q B*, 
B*B q (dashed) and B*B* (q = d,s), all generated with the same relative rate. In writing 

BqBq the first B q *^ is intended as the reconstructed one. 



one, due to reconstructed B mesons originating directly from the T(5S). This is 
shown in FigUl where the two components are independently shown in the ttles and 
m miss projections. 

We consider only VV events and we study the cross-feed generated by the tails of 
similar events (in this case Bd iS J/ip4>) coming from the other configurations (PP and 
PV). We have chosen specific selection criteria for ities or m m i SS to separate the different 
components. We find that in general, and mainly for the channels with photons in the 
final state as shown in the following sections, the m m i SS variables provide better efficiency 
and purity. In Tab. [2] we summarize the cuts applied on m miss , the efficiency and we 
quantify the cross-feed among B s and B d samples. 

The situation is different in the case of BBtt continuum production. These events are 
overlapped to B s mesons from VV events, in such a way that a simple cut in the AE vs. 
rugs plane cannot separate the two components without affecting the distribution of qq 
background. 

The separation is indeed possible since these events are characterized by a continuum- 
like distribution in tues (see left plot of Fig. [5]), while they exhibit a peaking structure in 
AE (see right plot of Fig. EJ), shifted to negative values with respect to B s mesons, because 



10 







5d 




B s 






S -> 






Cut 
Efficiency 


5.355 < 
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89% 
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98% 
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1.6% 
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B -> K*(K° s n°)j 


Cut 
Efficiency 


5.355 < 


m^ ss [GeV] < 5.410 
75% 


5.440 < 


m s miss [GeV] < 5.490 
84% 
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4.4% 




0.5% 




[PV= 


=3.3%, PP=1.0%] 


[PV= 


=0.4%, PP=0.1%] 



Table 2: Cuts, efficiencies and contaminations in the channels B —>■ J/ip<p ( no photons), 
B — > Kgir (2 photons), B — ► K*7 photons). The quoted values take into account the 
different production rates of the T(5S) resonance. The contamination of Bd to B s (B s to 
Bd) is negligible (less than 0.1% j. 

of the undetected energy of the pion. At present, concerning the relative amount of this 
component with respect to B s mesons, we can rely only the upper limit (UL) quoted by 
CLEO [15] (see Tab. [Q, so that we are not allowed to neglect this kind of background. 
Anyhow, the differences in the kinematic variables allow to isolate these events from the 
signal, adding a component in the Maximum Likelihood (ML) fit as it is currently done 
with BB background in many charmless analyses in BaBar [22]. 

3.2 Selection of B d s — ► K° s tt° 

As mentioned in Sec. |3j the presence of neutral particles in the final state smears the 
distribution of the interesting kinematic variables. The ability of the cuts to separate the 
different channels is therefore reduced by this effect, so that the result of the previous 
study is not guaranteed to hold in this case. To investigate this point, we consider a 
sample of B d , s K^tt decays, reconstructing signal candidates from — > and 
7T° — > 77 decays. 

Fig. [6] shows the niEs, remiss an d m miss distributions in this channel. Although the 
width of the distribution increases in a significant way, it is still possible to separate 
the different contributions. In particular, the use of m^ iss rather than m E s provides a 
separation almost as good as in the case of J/iftcf) final states. 

3.3 Selection of B d s — ► K*^ 

As a further step, we repeat the previous study using Bd, s — ► K*(Kg7i )'y signal 
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events, reconstructed from K s — > 7r + 7r~ and 7r° — > 77. In this case, there are three 
photons in the final state. As in the previous section, an improvement in the separation 
between different channels can be achieved using m miss instead of rriEs- In Tab. [2], we 
summarize the cuts applied on m miss in the three channels, with efficiencies and the 
details on the contaminations. 



3.4 A Remark on B reconstruction 

We would like to highlight a point that is relevant for the rest of this paper. Since 
it is possible to isolate B s mesons from and charged B mesons, it is clear that it is 
possible to continue the T(4S) physics program when running at the T(5S). This is 
particularly crucial for measurements that are difficult to perform at other facilities, such 
as inclusive measurements (b — > S7, b — > sll, and semileptonic decays) as well as some 
exclusive measurements with open kinematic (B — ► tv, exclusive semileptonic decays). 
All these measurements represent important milestones of the physics program of the 
current B- Factories and it is hard to imagine a precision test of the flavour sector of NP 
without the possibility of reducing the errors on these observables with future facilities. 
At the same time, as already stressed in Sec. [2j one should keep in mind that there is a 
price to pay in terms of statistics, corresponding to a factor about three (six) for inclusive 
(exclusive) measurements. 

4 Accessing CP asymmetries in B^—B^ events 

In order to cover the physics program of a i?-Factory running at the T(45) resonance, 
not only one should be able to identify and charged B mesons, but also to perform 
time-dependent measurements of CP asymmetries in B^ decays, which are based on the 
fact that the initial BB pair is produced in a coherent state. Moving to the T(5S) does 
not bring problems in terms of vertexing and tagging, since the experimental environment 
is similar to that of BaBar and Belle. In addition, as for the T(4S), the B^B^ pair is in 
a coherent state when the T(5S) decays in a PP final state. We show below that this is 
also true for the dominant contribution, namely VV decays. We also discuss how, using 
time-integrated measurements, PV decays could be used to access the value of Tm\ CP , 
with \ f CP = 

This last feature represent a novelty with respect to the study of T(4S f ) decays and it 
allows to obtain an information equivalent to the measurement of the coefficient S of the 
time-dependent CP asymmetry even for those decays for which a measurement of the 
B decay vertex, needed for time- dependent measurements, cannot be performed. This 
is for instance the case of B^ — > 7r°7r° and B^ — > 77 decays, for which otherwise only 7 
conversion processes would allow to try a time- dependent measurement, paying a big 
price in terms of reconstruction efficiency. 

4.1 The coherent time evolution of the B*B* mesons 

As explained in Sec. [2j at the T(5S) resonance B*B* meson pairs are produced in a 
J CP = 1 state. The question is whether the B q B q pairs resulting from the B* — > B q ^ 
and B* — > B^pi decays exhibit the same properties of quantum coherence of the initial 
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state. Since the B* — > B q ^ decay (and its CP conjugate) is an electromagnetic decay 
and can be considered instantaneous, one is allowed to regard the B q B q / -f-f final state 
as a direct product of the T(5S f ) resonance. This implies that the state constituted by 
these four final particles has to preserve the initial quantum numbers. Since the two 
photons have to satisfy the Bose-Einstein statistic, their sum has to be a symmetric state. 
This request forces the B g B g pairs to be in an antisymmetric state. This implies that 
B q B q pairs produced in a VV event of a T(5S f ) decay have the same quantum coherence 
properties of B°B° pairs generated at the T(4S) resonance. The full expression for the 
decay amplitude is: 

A» = B a ( Pl + k 1 )'B 8 (p 2 + k 2 )-G 3 - 

/e(q,p 1 ,k 1 ,e 1 )e(q,p 2 ,k 2 ,e 2 )p^ e(q,pi, fci, e 1 )e(q,p 2 , fo, e 2 )p 2 . 



2M 2 B , 2M 2 B , 



6(g,Pi,fci,ei)e(g,p 2 ,fc2,€2)fci _ e{q, pi, fci, ei)e(g, p 2 , k 2 , e 2 )k^ - 
2M%, 2M 2 * 



+ B s (p 1 + h) ■ B s (p 2 + k 2 ) ■ G 2 - 

■ ( - ^Q,Pi,ki,e 1 )e(p 2 ,k 2 ,e 2 ,fi) + e(q,p 2 , k 2 , e 2 )e(pi, fci, e u fi)+ 
+ 2e(pi, k 1 ,e 1 , cr)e(p 2 , k 2 , e 2 , a)rf - 2e(pi, fci, e u <r)e(p 2 , k 2 , e 2 , a)p%+ 



where 



+ 2e(pi, fci, ei, cr)e(p 2 , k 2 , e 2 , cr)fcf - 2e(pi, fci, ei, (r)e(p 2 , fc 2 , e 2 , cr)fc^ + 

+ B s (pi + fc 2 ) • 5 s (p 2 + fci) • G 3 - 
-e(q, p u k 2 , e 2 )e(q, p 2 ,k u e^pf e(g, p 1? fc 2 , e 2 )e(g, p 2 , fci, ei)^ 



2M B% 2M BJ 

_ e(g,Pi,fc2,€2)e(g,P2,fci,€i)fci e(g,pi, fc 2 ,e 2 )e(g,p 2 ,fc 1 ,ei)A: 2 ' \ 
2ilf2 + 2M2 y + 

+ 5 s (pi + fc 2 ) • B s (p 2 + fci) • G T 

■ ( - fa, e2)e(p 2 , fci, ei, //) + e(g,p 2 , fci, ei)e(pi, fc 2 , e 2 , /i) + 

+ 2e(pi, fc 2 , e 2 , cr)e(p 2 , fci, ei, cr)^ - 2e(pi, fc 2 , e 2 , cr)e(p 2 , fci, ei, cr)^- 

- 2e(pi, fc 2 , e 2 , cr)e(p 2 , fci, ei, cr)fcf + 2e(pi, fc 2 , e 2 , cr)e(p 2 , fci, ei, cr)fc 2 ) , ( 2 ) 



B s( P ) P 2_ M 2 +iM F 

e(a, 6, c, d) = ea^c^b 13 c' d s = € al3lS a a bpCjds 

e(a, 6, c, cr)e(d, e, /, a) = e a/3ia e SXuo 'a a b^ c° 'd 5 e x fu, (3) 

and /x is the index contracted by the T(5S) polarization vector, q is the T(5S) momen- 
tum, pi(p 2 ) is the B S (B S ) momentum, fci and fc 2 (ei and e 2 ) are the photons momenta 
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Figure 7: Statistical error on the measurement of Acp(Bd — > /or PV decays collected 

at the T(5S), as a function of the integrated luminosity. 



(polarization vectors), and G2 and G3 are two hadronic form factors [23] q Among 
the variables entering Eq. (j2J) and Eq. (j3J), G 2 and G 3 are the only whose knowledge is 
not precise. Different estimates exist in literature [23,24], but at present the theoretical 
expectation is controversial. Anyhow, since the ratio of the two determines the angular 
distribution of the B*B* mesons (and so also of the B q B q pairs) in the T(5S) rest 
frame, the study of the angular distribution of the final states in T(5S) decays allows to 
experimentally ascertain them. In any case, whatever the value of the two form factors 
is, the coherence of the B q B q is guaranteed by Eq. (J2J). 

4.2 Time- integrated CP asymmetries hi B*B events 

When the T(5S) decays to B*B events, the final BB system is in a C= +1 state, 
after the decay B* Bj. This difference implies that, unlike the case of T(4S*), for each 
value of the time t, before one of the two B mesons decays, the two mesons are both Bh 
or Bl CP eigenstates. Integrating the time out of the wave function of the BB pair, the 
time-integrated CP asymmetry becomes [25]: 

A f _fl-y 2 \ 2 (l-x*)(l- |A£ P | 2 ) + 4xlm(\f CP ) 
CP V1 + *V (l + y 2 )(l + \\ f CP \ 2 )-4y1le(\ f CP ) 

where x = Am/T, y = Ar/2r, \ CP = p^J, q/p is the mixing parameter of the B—B 

mixing, Af (At) is the amplitude for B — > / (B — > /) decays, and T]f is the CP eigenvalue 
of the final state /. 

While this might not be relevant for those channels for which X^ CP can be determined 
at the T(4S*), it opens new perspectives in the study of CP asymmetries for Bd decays 

3 The Gi form factors are defined from Eq. (3.2) of Ref. [23] after imposing the relation G\ = -2 G2, 
which is valid for neutral mesons in the limit where the Q 2 dependence in the charge form factor can be 
neglected. 
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Figure 8: Determination of the angle a without (top) or with (bottom) the constrain 
from direct CP asymmetries of Bd — > titt decays in PV events, combining 1 ab~ 1 (left), 5 
ab^ 1 (center) and 30 ab^ 1 (right) integrated at the T(5S) for 2 ab^ 1 collected at the T(4S) 
(end of BaBar and Belle data taking). 



to neutral particles, among which two remarkable examples are B^ 77 and B^ — > 
7T°7r . B q — > 77 decays are sensitive to NP effects in b — > q transitions and they are 
complementary to the measurement of b —>■ q'j processes (see Sec. I5.3.3[) . We do not 
discuss this specific case, while we concentrate our attention to the impact of measuring 
the direct CP asymmetry of B^ — > 7r°7r° decays in PV events. 

Measurements of rate and asymmetry of B^ — > 7r 7r are currently used in the isospin 
analysis of B — ► ttti decays for extracting a. The isospin analysis provides eight solutions, 
corresponding to the trigonometric ambiguities of the isospin construction [26] . Recently, 
it was shown that the use of basic QCD properties allows to exclude some of the solu- 
tions [27], still leaving several ambiguities in the determination of a. The measurement of 
Im(App) from Eq. (j3J) allows to remove some of these ambiguities, providing an additional 
information with respect to BR's and CP asymmetries of the other channels. 

From an experimental point of view, the determination of the direct CP asymmetry 
in PV decays is equivalent to what is done at the T(4S) for the same decay channel [28]. 
Using kinematic cuts and requirements on the energy deposits in the electromagnetic 
calorimeter, signal events are isolated from continuum background in a ML fit that uses 
topological and kinematic variables. Using the information obtained applying the tagging 
algorithm on the rest of the event, the flavour of the decaying B meson is determined 
on a statistical basis, taking into account the misidentification probability of the tagging 
algorithm. Scaling the currently available BaBar result to the statistics available at a 
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B-Factory machine running at the T(5S) and assuming the same tagging performances, 
we obtain the experimental error as a function of the luminosity shown in Fig. [7J 

To illustrate the impact of this measurement on the isospin analysis, we implemented 
this new measurement in both — > 7t 7t° and — > 7r + 7r _ modes in the analysis by 
UT fit [27]. The measurements currently available will be improved using VV (PV and 
VV) events for CP asymmetries (decay rates). We assume also that the results of the 
T(55) runs will be combined to the outcome of 2 ab _1 collected by BaBar and Belle at 
the end of their data taking. At the same time, the measurement of direct CP asymmetry 
for neutral B decays in PV events will provide the experimental determination of the 
observable in Eq. (J4]). While Acp(7r°7r°) provides a new information, with respect to 
what is currently used in the isospin analysis, Acp{ji + ^) contributes to reduce the error 
on Zm(Acp)(7r + 7r~), already constrained by S(tt + 7t~). For all the measurements, we 
use the current central values and the errors are obtained scaling the statistical errors 
to the assumed luminosity, without reducing the systematic error. The error on the 
measurements of direct CP asymmetries in PV decays are obtained scaling those of PP 
decays to the ratio of relative production fractions (see Tab. [I]). 

The result in terms of the determination of the angle a is shown in Fig. [8] using the 
current isospin analysis (top) and adding the information from PV decays (bottom), 
for three different values of integrated luminosities. The plots clearly show that the 
ambiguity is substantially reduced already with 5 ab" 1 of integrated data. 

5 Accessing the B®—B® mixing phase 

One of the main ingredients of the success of the current 5-Factories is the possibility 
of using the coherence of the initial B B state to exploit the interference between CP 
violation in the decay and in the mixing. Thanks to this technique the measurement of 
sin 2f3 has been possible with a high level of accuracy. 

In the same way, one would like to exploit the coherence of the B*-B* final state at 
T(5S). The paradigmatic physics goal in this field is the determination of the weak phase 
(3 S using B s — > J /ip4> decays. 

As described in the following, this does not seem achievable with the precision of the 
current vertexing detectors and the design parameters of the next-generation 5-Factories, 
currently under discussion [13]. The main problem comes from the large value of the ratio 
Am s / Airirf = 35.0 ± 0.4, which requires a presently unreachable vertexing resolution to 
be sensitive to B s oscillations. 

Nevertheless, information on the mixing phases can still be obtained using time- 
integrated measurements. In general, comparing the tagged decay rates for B s decays 
into CP eigenstates / for positive and negative values of At, it is possible to determine 
the value of TZeX^p and ImX^ p . Considering different decay modes, several combinations 
of weak phases can be extracted. In particular, we consider the case of the determination 
of p s from tree- level B° — > J/ip<j) decays and from penguin dominated B s — > K°K° 
decays. In addition, complementary information on f3 s can be obtained from the charge 
asymmetry in flavour specific final states or in dimuon events, or from the angular 
analysis of B° — > J/ip(f) decays. 
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5.1 Time-dependent CP asymmetry 

The study of time-dependent CP asymmetries is one of the milestones of the physics 
program for currently running B- Factories. It is based on the simple idea that, provided 
that the laboratory frame is sufficiently boosted with respect to the CM frame of the T 
resonance, it is possible to access the B q oscillation even though the vertex resolution is 
smaller than Am q , thanks to the (3^ Lorentz factor [29]. In the case of Bd-Bd oscillation, 
a boost of /?7 ~ 0.5 was large enough to allow the measurement of sin 2(3 with the vertex 
precision available when BaBar and Belle were designed. With the currently available 
technology on silicon vertex detectors, the boost could be in principle reduced by a factor 
of two or the resolution improved by the same factor using the current boost. Nevertheless, 
this is not enough for B s -B s oscillations, which are too much faster than in the case of 
Bd mesons. 

We used Monte Carlo simulated experiments (toy Monte Carlo experiments) to identify 
the minimal vertexing resolution needed to detect B s -B s oscillations at a 5-Factory. We 
quantify the minimal resolution as the largest value of the At resolution that allows to 
extract the S and C parameters of the time- dependent CP asymmetry of B s — > J/ip(f) 
decays with an accuracy comparable to the precision on sin 2/5. In this way, it is possible 
to provide an answer to the problem without necessarily relying on a particular set of 
machine and detector parameters. The result obtained can then be translated into a 
requirement on a given machine, once the energies of the beams (that fix the Lorentz 
boost) are defined. In fact, the precision on a (At) can be written in terms of the resolution 
of the vertex detector o~(Az) and the Lorentz boost /?7 of the CM frame with respect to 
the laboratory frame, by using the relation a(Az) ~ a(At)f3^c. 

We determine the shape of signal and background using fully simulated samples, se- 
lected requiring [30]: 

• 3.06 GeV < mjj^ < 3.14 GeV for J/ip -> //+//" 

• 2.95 GeV < mj^ < 3.14 GeV for J/ip -> e+e" 

• 1.004 GeV <m <s> < 1.034 GeV 

• particle identification (PID) requirements based on Cherenkov angle and dE/dx, 
both on J/tp and <j> daughters. 

• a cut on the polar angle of the thrust axis: | cos(0Thr) I < 0.8 

in order to reduce contamination from qq background and from tracks combinatoric. We 
also apply the following kinematic requirements: 

• 5.440 GeV < m s miss < 5.490 GeV 

• -0.2 GeV < AE < 0.1 GeV 

to select an almost pure sample of T(5S) — > B*B* events. This selection has a global 
efficiency of ~ 14.7% on signal events. We generate and fit samples containing signal 
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and background events (scaled to the considered luminosity), according to the following 
likelihood function: 



„ I , — r-r ' { ^sig ' Psig(jVlmiss)i ' Psig^AE)^ ' Psigi^-F)i ' -Psig(At)i 

y{N sig + N bkg )\ . =1 

Nbkg ■ Pbkg{ P bkg {AE) t ■ P bkg {T)i ■ P bkg {At)i}, (5) 

where P S i g (x) (P bkg (x)) is the probability density function (PDF) of the variable x for 
the signal (qq background) component. Background events from other BB decays are 
negligible after the selection cuts. P sig (At) takes into account the knowledge of the proper 
time of B mesons, coming from the vertex reconstruction. The full expression uses the 
convolution of the angular distribution function with the resolution function associated 
to the vertex reconstruction (usually, the sum of three Gaussians, with mean and width 
scaled according to the per-event error a At)- We assumed the models adopted for BaBar 
time-dependent measurements [31], for both the shape of At and the performances of the 
tagging algorithm. 

All the other functions are used to separate signal from background. Their shapes 
are obtained from unbinned maximum-likelihood fits to signal and background samples, 
coming from full Monte Carlo simulations. 

The signal distributions for m s miss) AE, and the Fisher discriminant T are parameter- 
ized as: 

/(x)=eXP [' 24Wx-m) 2 ] (6) 

where m is the maximum of the distribution, while a± and a± quantify the width and 
the tail, the positive (negative) sign corresponding to positive (negative) values of x — m. 
The invariant masses of the two resonances, and mj/^ are described by relativistic 
Breit-Wigner functions. 

In order to describe the qq background, we use a second-order polynomial for AE, 
while T is given by the function of Eq. (jH]) and m s miss is parameterized by a phase-space 
threshold function [34]: 



-(N sig +N bkg )/N x 



f(x) =aVl -x 2 exp ^-f • (1 -x 2 )J ; x = m s miss /m T (7) 

where tut is the value of the threshold and £ is the parameter determining the shape. For 
At we use a Gaussian resolution function. 

The Bd background is parameterized using the function of Eq. ([6]) for AE, and the 
phase-space threshold function (as in Eq. <^) for m s miss . All the other PDF's are similar 
to signal ones. 

For simplicity, we assume a fully polarized final state, which allows us to discard 
from the likelihood the PDF describing the angular distribution, and we parameterize 
the resolution function with a single Gaussian function, centered around —0.2 ps _1 (to 
include the effect of the charm bias on the tag side). Here, we do not use the per-event 
error but an average value a At as the RMS of the resolution function, varied in the range 
[0.02,0.20] ps _1 . We assume an integrated luminosity of 30 ah~ x and we generate and fit a 
set of toy Monte Carlo samples for each value of a At- The result of this study is reported 
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Figure 9: Error on S(J/if}<f>) as a function of the At resolution a^ t - Error bars represent 
the RMS of the error distribution for a given set of toys. 



in Fig. for the error on S. The error on C is not shown, since the values of the errors 
on S and C are the same for all the chosen values of oa*- 

It is clear from the plot that, in order to achieve an acceptable precision on the value 
of the CP parameters, values of cr(At) < 0.11 ps are needed. We also found that for values 
larger than 0.11 ps not only the average error increases, but also a bias is introduced in 
the fit, generating a larger spread in the error distribution (the error bars in the plot). 
For these values of a (At) the B s oscillation is too fast to be detected. The bad news is 
that the values of resolution required are beyond the possibilities of the current designs. 
The existing I?-Factories are outside the range of the plot, having a(At) ~ 0.65 ps. The 
improvements coming from new technology, together with the possibility of adding to 
the vertex detector a layer [8] on top of the beam pipe, allows to push a(At) down to 
~ 0.4 ps for a boost of ~ 0.28 [35], corresponding to ~ 0.2 ps for a boost of fij ~ 0.6. 

We stress the fact that statistics alone does not allow to reduce the errors in a 
significative way, since the limiting factor is in the vertexing resolution and does not have 
a statistical nature (as for instance in the case of decay modes with large backgrounds). 



5.2 Tagged rates at positive and negative At 

Even assuming that the needed time resolution for the study of B s -B s oscillations is 
outside the capability of the next-generation I?-Factories, it remains possible to measure 
weak phases in B s decays measuring the sign of At for B s -B s meson pairs. In fact, in this 
case one can measure decay rates as a function of B s flavour and At sign, which provide 
four experimental observables, depending on the weak phase of the chosen final state. 
Using this technique, it is possible to access the B s -B s mixing phase (3 S in B s — ► J/ijj(f) 
decays, as well as combinations of this phase and the phase 7 of the CKM matrix (as for 
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the case of B s — > D S K, which we do not discuss here). In this section, we introduce the 
basic formalism and we provide two examples of such measurements. 

5.2.1 Theoretical Formalism and Experimental Details 

Let us consider a B s pair produced at the T(5S) resonance, through a B*B* state. If 
one of the two B s decays into a CP eigenstate with eigenvalue rjf and the other one in a 
flavour tagging final state, the PDF of the proper time difference At can be written as: 



-|At 

P(At) oc e~ 



, fAT s At\ . . . , /Ar s At\ 

Ki cosh I — - — 1+^2 cos ( Am s At) + k 3 smh I — - — J + k 4 sin ( Am s At) 



where the coefficients are defined as: 



n, = -Ke\f CP ■ n 4 = - qtag lm\l CP . (9) 



These four quantities are effective parameters that depend on the tag sign q tag (+1 (—1) 
when the tag meson is a B s (B s )) and on the CP parameter \^ P . We use, as a good 
approximation, q/p = e~ 2l/3s in the following. 

In the case of the current B- Factories, it is usually imposed that Ar ~ 0, which 
holds with a very good accuracy for B^ mesons. This assumption reduces the number of 
parameters to two (the usual S and C). In the case of time-integrated measurements, 
the tagged rates are measured without any requirement on At, which is equivalent to 
integrating Eq. (jSJ) between — oo and oo. In the case of time-dependent measurements, 
the convolution between Eq. (jSj) and the resolution function of the vertex detector is used 
in the fit, which allows to extract the S and C coefficients. 

As discussed in the previous section, it is hard to imagine that time- dependent mea- 
surements of B s decays will be performed at a B- Factory. On the other hand, the presence 
of At-odd terms introduces an asymmetry between the numbers of events with At > and 
At < 0, depending on a£. p and q tag - In addition, the sign of At can be determined with 
good precision even without improving the vertex resolution, since tb s ~ T Pd . It is then 
possible to extract the four parameters from the tagged decay rates for positive and 
negative values of At, the dependence on the being provided by the integral of Eq. (|S]) 
and the resolution function of the vertex detector in the range [— oo, 0] and [0, oo]. As an 
alternative, one can measure a flavour integrated rate and quote flavour asymmetries (as 
currently done by the 5-Factories) , which allows to cancel some of the systematics. 

Exploiting the technique described above, it is possible to measure the parameter \^ CP 
for all CP eigenstates of the B s , without any need of improving the vertex resolution. Al- 
though this procedure is not as powerful as the time-dependent measurements performed 
by experiments at hadron collider, it can be applied to all the channels for which the At 
sign can be detected. Given the very clean environment of e + -e~ machines, the set of 
accessible channels is larger than in the case of LHCb. In particular, since the position 
of the B decay vertex can be measured even using Kg — > 7t + tt~ rather than prompt 
tracks [21,36], the possibility of accessing A^ p even in these cases opens up interesting 
possibilities. 
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Figure 10: Distribution of the error on (3 S from B s — > J/ip<f) fl e ft) an d B s — > K°K° 
(right) decays, as a function of the integrated luminosity, obtained from the measurement 
of tagged rates with At > and At < 0. For the right plot, the estimate of the theoretical 
error is also given (see text for details). 



Covering all the possibilities goes beyond the purpose of this paper. Nevertheless, we 
give below two examples of the potentiality of this new technique. 



5.2.2 Determination of (3 S from tree-level processes 

The study of B s decays at hadron colliders allows to determine the absolute value and 
phase of the B s -B s mixing amplitude. The comparison of the experimental measurements 
to the theoretical expectations allows to test the presence of NP in b — > s transitions. 
The recent measurement of Am s [10] already provided the first milestone of this physics 
program. A first measurement of the mixing phase is available from D0 [37], obtained 
from a time-dependent angular analysis of B s — > J/ip(f> decays. 

A large improvement will come from LHCb [38]. In fact, this experiment is expected 
to measure, with the same technique, both the width difference and the phase up to 
cr(Ar s /r s ) = 0.0092 and cr(f3 s ) = 0.023 [39] in just one year of nominal data taking. 
With the full expected B physics dataset (10 fb^ 1 ) the error on (3 S can be reduced to 
cr(/3 s ) ~ 0.01. It is clear that such a precision cannot be achieved with the Lorentz boost 
available at present e + e~ facilities and the resolution of actual vertexing detectors (see 
Sec. 15.11) . Nonetheless, we will show how a 5-Factory running at T(55) could provide 
alternative and independent informations on the same quantities, even if with not the 
same accuracy. 

In fact, one can use the approach described in Sec. 15.2.11 to measure the four 
parameters in the case of B s — > J/ipip and then extract from them a value for (3 S . As for 
the measurement of (3 from Bd — > J /if) Kg, one can neglect CKM suppressed subleading 
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contributions □ and assume that X f CP = Vf e^. 

In order to evaluate the sensitivity to (3 S , we used the same framework described in 
Sec. 15.11 to perform a set of toy Monte Carlo experiments, as function of the integrated 
luminosity. The statistical error as a function of the luminosity, shown in the left 
plot of Fig. [101 clearly proves the potentiality of this technique and the possibility to 
measure f3 s with a relatively good precision. We found a two-fold ambiguity between 
(3 S and —/3 S . When the value of . /3 S is close to zero (as it should be in the SM), two 
different and partially superimposed peaks appear, which produce a (almost) two-times 
larger resolution in the total pdf. Anyway, it is important to stress that a similar ap- 
proach can be followed for B s decays into higher charmonium resonances, which provides 
independent determinations of (3 a , taking into account the different hadronic uncertainties. 

5.2.3 Determination of f3 s from penguin modes 

The same experimental technique can be applied also to determine j3 s from penguin 
modes. The comparison of this result to the value measured in B s — > J / ip<f> allows to test 
the consistency of the SM and to constrain NP models. 

A similar test is usually performed in the Bd sector comparing the value of sin 2f3 from 
Bd — * J/i/)Kg to the S value of penguin dominated modes, such as B d — > (j)Kg. The 
equality {sa\20)j/^ K o = (sm2j3) penguin would be strictly true only if these decays were 
mediated by a single combination of CKM matrix elements, while in all the cases a CKM 
suppressed amplitude is present. At the very high precision that is expected with the next 
generation of B physics experiments, it is not be possible to neglect the CKM suppressed 
contribution, which introduces a theoretical error associated to the SM expectation of 
AS = S(b^ sqq) -sin 2/?. 

The same problem is present in the case of B s — > K°K°, since the decay amplitude in 
terms of renormalization group invariant parameters is given by the relation [42]: 

A(B S -> K°K°) = - V us V: b P? 1M -V ts V t * b P s , (10) 

which is formally equivalent to the amplitude of the golden mode B^ — > <f)K°. The main 
advantage of considering this decay is that the theoretical error can be estimated using a 
data-driven approach [40,41], as we explain below. 

As a first step, we assume that P S GIM can be neglected, which implies that \(, P = 
rjjre~ 2l ^ s . We use a set of toy Monte Carlo experiments to extract the expected error on 
the weak phase, applying the method described in Sec. 15.2.11 We determine the error on 
(3 S as a function of the integrated luminosity, as shown in the right plot of Fig. [TD1 

Comparing the two plots of Fig. [101 it is clear that B s — > K°K° is affected by a larger 
background, which increases the error on the weak phase, the dominant contamination 
coming from qq events (q = u,d,s,c). These events are characterized by a jet-like dis- 
tribution in the center of mass of the e + e~ system, which allows to separate them from 
signal events. At a new facility, the discriminating power might benefit from the improve- 
ment of the vertexing resolution. Using modern technology, it is possible to improve the 
vertexing precision [43] , allowing to separate the B and D vertices on the tag side of a BB 
event. Since no secondary vertex is present in qq events, it is not unrealistic to imagine 

4 As it is done for Bd — > J/ipK®, [40] , the impact of the CKM suppressed subleading contributions can 
be obtained from the measurements of A^ P for similar channels, such B s — ► J/tpK*° 
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that such an approach will strongly suppress the background contamination, reducing the 
error on f3 s . A quantification of this improvement relies on a specific detector design and 
goes beyond the purpose of this paper. 

The estimate of the error induced by neglecting P s GIM can be obtained considering the 
time-dependent study of B d —>■ K°K°. The decay amplitude is given by 

A(B d ^K°K°) = - V nd V: b Pr-V td V; b P d , (11) 

which is similar to Eq. fflOl) . except that here the two combinations of CKM elements 
have the same order of magnitude. This difference maximizes the sensitivity to the ratio 
P d IM /P d , which can be determined assuming the SM values of the CKM parameters and 
using the experimental values of the BR and the CP parameters S and C to constraint 
the hadronic parameters. These hadronic parameters cannot be assumed to be equal 
to those of Eq. (fTOj) . because of SU(3) breaking effects. Nevertheless, it is true that an 
SU(3) breaking larger than 100% has never been observed up to now. Thanks to this 
consideration, one can determine the maximum allowed value of P^ IM /P s , taking for the 
output distribution of P d IM ' / P d (to take into account the statistical error on the fit to 
Bd — > K°K°) a 100% range centered around the mean value of P^ IM jP d (to take into 
account SU(3) breaking effects). From this determination, an estimate of the induced 
error on f3 s can be obtained, in analogy of what is done for /3 from Bd — > J/i/jK° using 
B d -»■ J/V>7r° [40]. 

In order to give an estimate of the theoretical error induced on (3 S by neglect- 
ing P{f IM /Pd, we scaled the statistical error of the time-dependent measurement of 
Bd — > K°K° by BaBar [44], assuming an irreducible systematic error of 0.01 on S and 
C . The output is shown in the plot on the right of Fig. ITOl as a function of the integrated 
luminosity. Here, we assume that data will be collected at the T(5S) resonance, but it 
is important to recall that, for a certain amount of integrated luminosity, the number of 
Bd decays available from T(4S) is larger (so that the error corresponding to the same 
luminosity will be smaller, if not systematics dominated). 

5.3 B s —B s mixing phase from time- integrated measurements 

In order to extract the weak phase (3 S of B s -B s mixing, one can use an independent 
strategy, which, as for the measurement discussed in Sec. 15.2.11 does not rely on 
time-dependent CP asymmetries. In fact, it is possible to use the determination of 
Ar s [47] and the charge asymmetry in semileptonic B s decays, A S SL [45] to obtain an 
independent determination of the same quantity. As we show in the following sections, a 
good precision on these two quantities can be obtained at a 5-Factory, exploiting high 
statistics, high efficiency in lepton reconstruction and a pure B s sample. 

5.3.1 Measurement of Ar s /r s 

The relative decay-time difference AT S /T S can be determined studying the angular 
distribution of the B s — > J / ip(p decay. Using the transversity basis, the angular distribu- 
tion of the final state can be described in terms of three complex amplitudes: the two 
transverse amplitudes A±_ and A» (perpendicular and parallel), and the longitudinal com- 
ponent A . A and A\\ are CP-even, while A±_ is CP-odd [46]. In the SM, neglecting CP 
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violation effects in the B s mixing, the CP eigenstates Bcp = +1 and Bqp = — 1 corre- 
spond to the mass eigenstates B L and B H . CP-even (odd) amplitudes evolve according to 
the exponential factor e~ Tljt (e~ FHt ). As suggested in ref. [47], if CP violation in mixing 
is not negligible (which is excluded in the SM but it is possible in NP scenarios) the time 
evolution of these amplitudes is modified. In particular, an explicit dependence on the 
CP violating weak phase (3 S appears. Following these considerations, one can write the 
angular distribution of B s events as: 



dpdt 
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^ 2 l + \cos(2p 8 )\ h ^ + |A |2 l + |cos(2&)| 



h (P) + 
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(12) 



where the functions /i(p) and the transversity variables p = {cos 9, <p, cos-?/'} are defined 
in [47], and <5i (<5 2 ) is the strong phase difference between A± and A\\ (A± and A ). The 
arbitrary phase in the decay amplitude is removed forcing A to be real. In the PDF, we 
write Th and Tl in terms of Ar = Tl — Th and T = (Tl + r#)/2, and we float these 
two parameters in the fit, along with the yields, the absolute values of the amplitudes, 
the strong phases S\ and <5 2 and the weak phase (3 S . In the MC generation we assume 
Ar s = 0.1 ps _1 (SM prediction from Ref. [32]), and Y s = 1.39 ps _1 [33]. It is important 
to stress the fact that no assumption on the size of CP violation in B s -B s mixing has 
been done, which is important to constrain physics beyond the SM. 

We determine the shape of the m miss and AE distributions in the B s — > J '/ip<f> sample 
by means of a full MC simulation, applying the same requirements of Sec. 15.11 Then, 
in order to determine the expected error on Ar s /T s and (3 S , we perform a set of toy 
Monte Carlo experiments, generating a large number of datasets (corresponding to the 
assumed luminosity, and taking the the value of reconstruction efficiency from simulation, 
and BR(B S — ► J/ip(f) = 9.3 ■ 10~ 4 ) and fitting for the value of Ar s /T s , together with 
signal and background yields, through an extended and unbinned ML fit, by using the 



5 Here Tl and Th represent the lifetimes of the light and heavy mass eigenstates of the B s -B s system. 
Notice also that at a B-Factory the time t can be substituted by At, the proper time difference between 
the decays of the two B mesons 
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likelihood function of Eq. (jSj), with the addition of the PDF for the angular distribution. 
For the signal component we use the expression of Eq. (1121) . The angular distribution of 
background events is assumed to be flat in p variables. 

We performed different sets of toy Monte Carlo experiments, corresponding to different 
integrated luminosities between lab^ 1 and 30a6 _1 . Fig. [TT1 shows the pull distributions 
for Ar s , T and f3 s at 30 a& _1 and the dependence of the statistical error on the luminosity. 

We remark here that the measurement can be affected by large correlations between 
the measured quantities, as experienced by the D0 Collaboration in a similar analysis [37]. 
Also these correlations can be studied in the toy Monte Carlo experiments, and they can 
be used to extract the corresponding constraints on the NP parameters (see Sec. |7J) 

5.3.2 Charge asymmetry in semileptonic decays 

The amplitude describing B s -B s mixing can be experimentally constrained looking at 
the difference between the B s and B s semileptonic decay rates after the mixing of the B s 
mesons. Since the charge of the lepton tags the flavour of the B s meson at the decay point, 
and since the decay amplitude is characterized by only one combination of CKM elements, 
any deviation of the charge asymmetry from zero can only come from a deviation of \q/p\ 
from one. One can define the semileptonic asymmetry as : 

A s _ BR(B S ^B s ^ PP-l+Uj) - BR(B S -> B s -> D^Ija) _ 1 - \g/ P \ 4 {u) 
SL BR(B S — > 5 S — > Dl*)~l + vi) + BR(B S —> B s ^ D^ + l~Vi) 1 + k/p| 4 ' 

To measure A S SL , we exclusively reconstruct one of the two B mesons into a self-tagging 
hadronic final state (such as B s — > Da n) and look for the signature of a semileptonic 
decay (high momentum lepton) in the rest of the event (ROE). Since both the B mesons 
are reconstructed from self-tagging modes, it is possible to isolate those events in which 
the two B s mesons have the same flavour content. To increase the statistics, one can also 
avoid to reconstruct one of the two B mesons and determine its flavour on a probabilistic 
basis, using a tagging algorithm. 

To evaluate the expected error associated to this strategy, we use a toy Monte Carlo 
technique, generating and fitting a large set of toy-simulated data samples. We use two 
variables to discriminate signal from background, namely the m s miss of the fully recon- 
structed B s — > D s *\ candidate, and the missing mass of the other B, m u . Signal events 
peak in the region m s miss ~ 5.46 GeV and m u ~ GeV. We describe the event distribu- 
tions using the function of Eq. (jSJ) for m s miss and the sum of four Gaussians for m v . The 
continuum background, coming from the combinatoric of particles in the hadronization 
of e + e~ — > qq events (q = u, d, s, c) (~ 10% of the signal yield) is described by the phase- 
space threshold function of Eq. (j7|) for rn s miss and a second order polynomial for m v . An 
additional source of background comes from other B s decays (~ 2.5% of the signal yield), 
which can be described by the m s miss PDF used for signal and an order one polynomial 
for m„. Background from T(5S) — > BBn events is strongly suppressed, since both the B 
mesons of the event are requested to decay into D s mesons, which occur through CKM 
suppressed amplitudes for B d and B + mesons. 

We calculated the expected qq background yield scaling the values obtained by 
BaBar [48]. For signal, we consider both D s lv and D*lv events (I = e,p), assuming 
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Figure 11: Distributions of the pull for AT S (top row), T s (middle row) and /3 S (bottom 
row) from a set of 900 toy Monte Carlo experiments, assuming an integrated luminosity of 
30 ab^ 1 (left) and the trend of the error as a function of the integrated luminosity (right). 
The error bars, in the right plots, show the RMS of the error distribution. 



26 



0.007 r-r- 

_l 

CO 

< 

^ 0.006^ 
0.005 7 
0.004 r 
0.003 r 
0.002 ~ 
0.001- 



10 15 20 25 30 

Integrated Luminosity[ab" 1 ] 



Figure 12: Statistical error on A S SL , as a function of the integrated luminosity. The 
dashed line represents the order of magnitude of the current systematic error on Acp 
measurements, shown for comparison. 



BR(B S -> DJu) = 2%, BR(B S -> D* s lu) = 4%, and the efficiency and purity of the 
BaBar analysis. 

The result of the simulation is shown in Fig. [121 where the statistical error on A$l is 
given as a function of the integrated luminosity. The dashed line on the plot represents 
the current value of the systematic error typically quoted in precision measurements of 
direct CP asymmetries at the 5-Factories. It is not unrealistic to expect an improvement 
of a factor two in the evaluation of the systematic error, considering that a sizable part 
of it has a statistical nature (being related to the available statistics for control samples). 



5.3.3 Inclusive Measurement of Dimuon Charge Asymmetry 

Running at the T(4f>), the measurement of A d SL can also be performed in an inclusive way, 
identifying BB events from pairs of same-sign leptons with positive (N ++ ) and negative 
(N ) charge and calculating the charge asymmetry as 

Ad - N+ + 

A SL ~ N — + N++ ■ ( 14 ) 

These events correspond to T(4S) decays in which both the B mesons decay into a 
semileptonic final state. Requiring the charge of the leptons to be the same, is equivalent 
to select only events for which the two mesons have the same b flavour content, which 
can happen only after a B B mixing process. Contrary to the analysis presented in the 
previous section (A S SL ), the cascade events component has to be controlled. A non-null 
value of A d SL is sensitive to CP violating effects in the mixing process, i.e. to the weak 
phase of the mixing. 

The situation of the T(5S) is more complicated: B s and B^ u mesons are produced 
simultaneously, but cannot be separated during the reconstruction, since the analysis is 
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done inclusively. In this case, the measured observable can be written as 



A 



N— -N++ x-X 



(15) 



CH 



N" +N++ x + X- 2XX 



where ^ = (1 - f s )xd + fsXs, X = (1 ~ fs)Xd + fsXs, fs is the production rate of B s meson 
pairs at the T(5S) and 



with z q = \q/p\ 2 q and z q = \p/q\ 2 q . 

The interesting feature of this observable is that it relates NP effects in and B s 
sectors, providing an additional constraint with respect to the measurement of Ag S L . The 
sample used for Aqh can be orthogonal to the one used for A S SL {A d SL ) measurements if 
one avoids the use of semileptonic decays to reconstruct the other B meson in the latter 
analysis. 

From the experimental point of view, the main background sources are qq events 
(q = u, d, s, c) and cascade-lepton events in which one of the two leptons is generated from 
the semileptonic decay of a D meson, coming from the B decay. These backgrounds can 
be separated from signal looking at the distribution of the decay time [49]. The charge 
of the muon distinguishes B and B mesons, while the information on the shape of the 
event in the CM frame allows to improve the rejection of qq background. The residual 
BB background, coming from cascade-lepton events can be suppressed using a Neural 
Network algorithm, based on: 

• the momenta of the two leptons with the highest momentum in the T(5S) center of 
mass system; 

• the total visible energy and the missing momentum of the event in the T(5S) center 
of mass system; 

• the opening angle between the leptons in the T(55) center of mass system; 

and trained on simulated samples of signal and background events. 

A rough estimate of the expected error achievable on the inclusive approach can be 
obtained scaling the statistical error quoted in the corresponding BaBar analysis [49] 
according to the inverse of Ng/\/Ns + Nb, where Ns and Nb are the expected yields for 
signal and background respectively. The result, as a function of the luminosity, is given 
in Fig. [13j The dashed line represents the order of magnitude of the current systematic 
error on Aqp measurement, shown for comparison. The measurement will be systematics 
dominated after a relatively small period of data taking, even though also in this case 
an improvement of the systematic error with an increased size of the control samples is 
possible. Notice that the information about A$l comes only from same sign events (i.e. 
with 2 leptons of the same sign), so only these events are taken into account in the error 
estimate. 
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Figure 13: Statistical error on Aqh, a function of the integrated luminosity. The dashed 
line represents the order of magnitude of systematic error on the current measurement at 
T(4S), shown for comparison. 



6 Benchmark Measurements of B s rare decays 

The production of B s mesons at the T(5S) allows to study the decay rates of the B s 
sector with the same completeness and accuracy that is currently available for the 
and charged B mesons, improving our understanding of B physics and helping to reduce 
the theoretical uncertainties related to NP-sensitive quantities. For instance, using the 
measurements of BR's of B s mesons to open charm and to charmless decays it is possible 
to close the SU(3) multiplets of b decays [50] and to control SU(3) breaking effects in 
channels like B& — > 4>Kg [51]. In addition, even without accessing time-dependent CP 
asymmetries, it is possible to access weak phases of B s decays. This is the case of the 
determination of 7 from the tree-level amplitudes of B s — > Kttit decays, for which it is 
possible to disentangle penguin and tree contributions using the distribution of the events 
on the Dalitz plot [52]. 

On the other hand, B s physics provides a direct way to probe NP effects in the b — ► s 
transitions. Effects generated from NP particles are expected to be more evident in b — *■ s 
than in b — » d transitions, not only because, thanks to the CKM hierarchy among the 
amplitudes, the effect of hadronic uncertainties on the calculation of the SM is smaller 
than in the case of b — > d processes, but also because this scenario is still allowed by the 
generalized UT analysis beyond the SM [32] and is expected in several extensions of the 
SM [53]. 

In the following we concentrate on this point. We show what could be learned from a 
high statistics study of B s mesons by giving some benchmark examples. In some cases, 
only at very high statistics and with particular detector performances, one can reach the 
same precision as at LHCb. Anyhow, one should take into account that the set of mea- 
surements we consider is not an exhaustive representation of the potential of a 5-Factory 
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running at the T(5S) and that a quantitative comparison between the performance of 
LHCb and a next-generation 5-Factory, based on a common set of measurements, goes 
beyond the purpose of this paper. What we rather want to prove is that, contrary to what 
is usually claimed, all the interesting measurements (rates of rare decays, angular anal- 
yses, measurement of direct CP asymmetry and weak phases in B s -B s oscillations) can 
be performed even if B s -B s oscillation cannot be accessed. To demonstrate this point, 
we already documented in Sec. [1] how the B s -B s mixing phase (3 S can be measured in 
tree-level and in penguin modes, testing the SM in B s -B s mixing as well as in b — > s tran- 
sitions. In this section, we discuss two other milestones of a B s physics program, namely 
the measurement of the BR of rare decays of the B s , and the test of the SM through the 
determination of |Kd/Ks| in penguin dominated modes. 

The main result of this study is that a S-Factory running at the T(5S) will not 
be limited by detector performances. Using the analyses techniques we briefly describe 
here, it will be possible to study several decays which might go beyond the capability 
of LHCb, such as channels without primary tracks in the decay and/or channels with 
only photons in the final state. Furthermore, the possibility of studying channels with 
open kinematic (such as b — > ulv and b — ► suu), even with a factor three reduction 
in statistics with respect to the T(4S) represents something that only at a B- Factory 
can be determined. In our opinion, the sum of all these arguments establishes the 
complementarity between the physics program of LHCb and the possibility of studying 
B mesons at the T(5S). So, while LHCb will provide high precision measurements 
for a limited set of observables, a B factory running at the T(5S) will be able to 
acccess (not nececssarely with a comparable precision) a much larger ensamble of 
quantities. If this is done before the LHCb will start collecting data, already with 
two ab _1 it will be possible to exclude NP contributions comparable to the SM 
amplitude, as well as to provide a set of measurements of rates and asymmetries to 
be used by LHCb as a calibration of the physics measurements in the very first run period. 



6.1 Probing New Physics in penguins with |Vtd/14s| 

The measurement of Am s at the Tevatron [10] has provided a further evidence of the 
consistency between experimental data and the predictions of the SM. The agreement 
between the prediction by UT fit [54] and the measured value of |K<i/Ks| has reduced the 
possibility of observing large NP effects in B s -B s oscillations [32] . Nevertheless, it is still 
interesting to probe NP in b — > s channels using an independent determination of |Vi<j/Vt s |, 
which might be sensitive to NP in a different way than Am s . It has been proposed in 
literature that such a test can be provided by the ratio R = BR(B° — > p 0/ y)/BR(B —> 
K* -f) [55]. This quantity can be expressed as 







V td ? 







where ttt4 is the combination of CKM elements one wants to determine: ^/i-^L, is 

\V ts \ ' |a£(K* u 7)|" 

the difference in the penguin factorized decay amplitude, due to the different mesons in 

2 2 

the final state; r m is defined as r m = ( ^ B __ ) 3 ; £ is the ratio of form factors F(B — > 

m B m K*0 

K*°)/F(B — ► p), and AR is a correction generated by the annihilation contribution to 



30 




Figure 14: Annihilation contribution to the b — > d amplitude of B — > p7 decays. 



B —>■ p7 (see Fig. [T41 . which is not present in B —>■ K* 0/ y. 

This ratio provides an independent test of the SM prediction of | Vtd/V ts \, since radiative 
decays are particularly sensitive to NP contributions. As these processes violate chirality, 
the decay amplitude in the SM is proportional to m q , the mass of the quark (g = b, s) 
generating the helicity flip. In presence of NP, the helicity flip is obtained by the mass 
insertion of a new heavy state, determining the enhancement of the NP amplitude with 
respect to the SM by a factor mHs/ m q, where itlhs is the mass of the new heavy state. 
This implies that it is possible to observe NP effects in these channels even when the 
contribution is too small to be determined in the B-B mixing. 

A limitation to the effectiveness of this test comes from the presence of the AR term in 
Eq. ffTT|) . In the SM this contribution is not only expected to be O^KQCD/^b) but it is also 
CKM suppressed, being proportional to cos(a) [55]. These two suppressions reduce the 
sensitivity to non perturbative QCD effects that determines the value of AR. The main 
limitation is related to the fact that, beyond the SM, the CKM factor is not necessarily 
small. While in MFV models the cos(a) suppression still reduces the sensitivity to the 
hadronic corrections, this is not the case for models with a generic flavour structure. This 
is why it is particularly interesting to look for a similar observable which is not affected 
by the presence of the annihilation term. 

Such observable can be provided by the ratio R = BR(B° S K*°^)/BR(B d -> K*°^y). 
In this case, in fact, the two decays are not affected by annihilation contributions and the 
ratio is governed by Eq. ffTTl) . with r§r^K-, and £ replaced by the similar factors, and 
where the AR term vanishes. 

The BR measurement for B^ — > K* 0/ y has been performed at the current B- 
Factories [56] and the present error is (almost) dominated by the systematic error. We 
assume BR(Bd — > K* 0/ ~f) = (40.1 ±2.4) ■ 1CT 6 , corresponding to the current world average, 
with the error given by the current systematic error. 

An estimate of the error that can be obtained for BR(B° S — > K* 0/ y) has been extracted 
performing toy Monte Carlo simulations as a function of the integrated luminosity. The 
SM expectation value for BR(B® — > K*°j) is about two times its B^ counterparts, £?° — ► 
p°7. As a reference, we use BR(B° S — > K* 0, ~f) = 1.8 ■ 10~ 6 , obtained using the central 
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Figure 15: Distribution of the statistical error on BR(B° S -> K*°j) from toy Monte Carlo 
experiments (left) and the expected error of the \V t d/V ts \ measurement (right), as a function 
of the integrated luminosity. 



value of the current world average for B® — > p°7 [57] and the relative factor two, coming 
from the Clebsch-Gordon coefficient of the p°. 

The distribution of the relative errors on the BR from the toy experiments as a function 
of the integrated luminosity is shown in the left plot of Fig. [151 The information coming 
from the toy experiments, combined to the current predictions on £ (£ = 1.17 ±0.09 [58]), 
allows to determined the value of 1^4. The result is shown in the right plot of Fig. [151 By 
repeating the exercise after reducing the error on £, we verified that the determination of 
is dominated by the statistical error even for large values of integrated luminosity. 

6.2 B s — > fji^L and the flavour structure of New Physics 

This decay plays a crucial role in study the structure of the Higgs sector and assessing 
the value of tan (3. In the SM, the BR of this process is written as [59] 

Br(B s -> i+r) = r(B s )^ ( " Y f% m?m B s \ 1 - 4^-\V£V ts \ 2 Y* (x t ) 7 (18) 

where x t = m^/m^Y, fg a is the B s meson decay constant and, within the SM, Y(x t ) is the 
Inami-Lim function [60], given by the expression 



Y(xt) = ? 



x t - 4 3x t ^ ^ 
_Xt - 1 (x t - l) 2 ' 



(19) 



The SM expectation for this decay rate is BR(B S — ► /i/i) = (3.35 ± 0.32) ■ 10~ 9 [61]. 
The decay occurs in the SM through loop diagrams, which make this process particularly 
sensitive to NP contributions. For instance, the observation of this decay above the 
SM prediction will provide a clear evidence of new heavy states contributing to virtual 
processes. Recently, a combined analysis of B and K rare decays [62] has addressed 
this point in a quantitative way in the context of MFV models at small tan /3, obtaining 
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BR(B S — > /i/i) < 7.4 • 10~ 9 at 95% probability. This result indicates that in MFV at small 
tan f3 we do not expect a large enhancement of the decay rate. Conversely, in very large 
tan f3 scenarios [63] , tan /3-enhanced terms can still give sizable contributions resulting in 
a enhancement of the decay rate well beyond the above limit. 

In order to estimate the expected experimental sensitivity in the SM, we assume the 
expectation value quoted above. We only consider B s -B s pairs from VV events. Since 
BR(Bd — > /i/i) is suppressed by more than a factor of ten with respect to BR(B S — > fifi), 
both in SM and in NP scenarios, we can neglect B^B^n events. The contamination 
coming from B s — > titt events is lower than in the case of B^, since the decay is CKM 
suppressed. As for the search of Bd — > fifi performed by the S-Factories [64], the hadronic 
background can be suppressed imposing a set of requirements on PID, typically combined 
into a PID neural net. The cost is about a factor 60% for the reconstruction efficiency 
With respect to the 5-Factories, we can also imagine a detector configuration with larger 
acceptance, thanks to a reduction of the Lorentz boost . Assuming a typical value of the 
reconstruction efficiency of 60% we expect 6 signal events and 960 background events in 
30 ab -1 of integrated data. 

We use a set of toy Monte Carlo to evaluate the sensitivity of this measurement using 
m miss , AE and the Fisher discriminant T as discriminating variables. We first generate 
signal and background events according to the shape of kinematic (m m i SS and AE) and 
topological (J 7 ) variables, as determined from fully simulated Monte Carlo samples. Signal 
events peak at a value of m m i SS ~ 5.46 GeV and AE ~ —47 MeV. Both distributions are 
described by Eq. (Q. For background events we use the threshold function of Eq. (J7J) for 
m miss and a second order polynomial for AE. For T, we use a bifurcated Gaussian to 
parameterize the signal distribution and the sum of two Gaussians for the background. 

For each generated datasets, we perform a ML fit as a function of the signal yield and 
obtain a shape of the likelihood for N^. We set the 90% upper limit to the value Nul 
such that 

f NuL C{N^)dN, 



90%. (20) 

)o ' w^*- 1 'MM 

We repeat this procedure for a large set of toys, as a function of the integrated luminosity. 
The result is shown in Fig. [161 where the full line represents the SM expectation and the 
errors bars the RMS of the upper limit distribution for each set of toy experiments. For 
comparison, the current 90% C.L. experimental UL (BR(B S — > /i + /i~) < 8.0 • 10~ 8 [65]) 
falls outside the range of the plot. 

Since in several NP scenarios we expect a sizable contribution to the decay amplitude 
from new heavy states, we repeated the exercise assuming a decay rate one order of 
magnitude larger than the SM prediction. In this case, it should be possible to measure 
the BR with a meaningful statistical significance, rather than setting an upper limit. 
This is shown in Fig. [THl where the statistical significance is given as a function of the 
integrated luminosity. 

As for the case of B s — > K°K°, discussed in Sec. I5.2.1[ this channel is affected by 
a large background contamination from qq continuum events (q = u,d,s,j), which are 
characterized by a jet-like structure, with all the particles coming from the primary 
vertex. Considering the possibility of an upgrade of the vertexing detector with modern 
technologies, it is possible to use the improved vertexing resolution [43] to separate the 
B and D vertexes on the tag side for signal events and reject qq events, for which no 
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Figure 16: Distribution of the upper limit at 90% C.L. (left) and statistical significance 
on the BR (right) for BR(B S — > fi + fi~), as a function of the integrated luminosity. Error 
bars represent the RMS of the UL distribution for each set of toys. For the left (right) 
plot, the SM expectation for the BR (ten times the SM value) was assumed. On the right 
plot, the dashed line represent the 5a evidence. 
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Figure 17: On the left upper limit distribution at 90% C.L. for BR(B S — > ^ ^r) in 
a Standard Model scenario, assuming a reduction of the background of a factor five as a 
function of the integrated luminosity. Error bars represent the RMS of the UL distribution 
for each set of toys. The horizontal band represents the Standard Model prediction. On 
the right significance of the measurement for BR(B S — ► in a New Physics scenario 

with an order of magnitude enhancement with respect to the Standard Model assuming a 
reduction of the background of a factor five, as a function of the integrated luminosity. 
Error bars represent the RMS of the significance distribution for each set of toys. The 
dashed line represent the 5a evidence. 



34 



secondary vertex is present. Using this additional requirement, it could be possible to 
strongly reduce the continuum background, improving the experimental precision. To 
give an idea of the impact of this improvement on the search of rare decays, we repeated 
the previous exercise assuming a reduction of continuum background by a factor of five, 
which might not represent an optimistic expectation^ We show in Fig. [17] the impact 
of this improvement on the result of our toy studies, for both the SM and the enhanced 
values of the BR considered above. In this case, one can obtain an observation of the BR 
in the enhanced scenario with a relatively small amount of statistics for a 10 36 cm _2 sec _1 
super 5-Factory. 

It is clear that in this case LHCb can perform much better, given the larger amount 
of B s produced in its environment and the clear signature of the muon tracks. In fact, 
LHCb is expected to be able to exclude the SM branching ratio already with 0.5 fb^ 1 
of collected data [39] and to reach 3a and 5a evidence of the SM signal respectively 
with 2 fb^ 1 (one year) and of integrated luminosity. Moreover, also ATLAS and CMS 
experiments are expected to give 4a signal after one year of nominal data taking at the 
luminosity of C = 10 34 cm _2 sec _1 . 



6.3 Measurement of B s — > 77 

For several years, b — > 37 has been considered as the golden mode to probe NP in the 
flavour sector. From an experimental point of view, any attempt to study this transition 
has to face the fact that the present error on BR(b — > 57) is dominated by systematic 
effects. It is then interesting to look for other channels that can play a similar role in 
constraining NP, but do not imply the experimental complexities that are introduced in 
b — > S7 by the knowledge of the photon spectrum. In this context, an interesting candidate 
is provided by B s — > 77 decays. The final state contains CP-odd and CP-even states, 
allowing to study CP violating effects with the tagging algorithms usually adopted at B- 
Factories. The SM expectation for the BR is BR(B S -> 77) = 2 - 8 • 1(T 7 [67]. NP effects 
are expected to give sizable contributions to the decay rate in some particular scenarios. 
For instance, in R-parity violating SUSY models, neutralino exchange can enhance the 
decay rate up to BR(B S — > 77) ~ 5 • 10~ 6 [68]. On the other hand, in R-parity conserving 
SUSY models BR(B S — > 77) is found to be highly correlated with BR{b — > 57) [69]. 
In fact, thanks to an extension of the Low theorem [70], the b —>■ 377 operator can be 
expanded at 0(Gf) on the standard operator basis needed for b — > 57. Doing so, one 
can relate deviations from SM in b — > 57 and B s — > 77 processes, provided the fact that 
NP effects appear only in the matching of Wilson coefficients of the standard basis of 
operators of the b — > 57 effective Hamiltonian. 

As already mentioned, the knowledge of the photon spectrum will always be a limiting 
factor for the measurement of b — > 37 decay rate, whose error is presently systematically 
dominated and it is not expected to improve in the near future. On the other hand, the 
exclusive measurement of B s — > 77 decay is very similar to other measurements already 
performed at 5-Factories (such as B° — > 7r°7r°) and it is not affected by a theoretical 

6 A better estimation of the impact of the vertex improvements on the data analysis can be obtained 
with a more accurate simulation of event reconstruction, which implies a specific design of the detector. 
This goes beyond the purpose of this paper, but it represents an interesting exercise to be considered by 
present and future experimental collaborations. 



35 



f 

i 

t 

t 



: 

5 10 15 20 25 30 

Integrated Luminosity[ab" 1 ] 

Figure 18: Statistical significance of the measurement of BR(B S — > 77) as a function of 
the integrated luminosity. Error bars represent the RMS of the significance distribution 
for each set of toys. 



error associated to the measurement (unlike the case of b — > 57 analyses). From an 
experimental point of view, the limiting systematic factor is represented by the knowledge 
of the efficiency for photon reconstruction, which can be reduced with dedicated studies 
on control samples with similar energy range for the photons. On the other hand, the 
presence of two hard photons represents a clear signature for signal events, in particular 
if studied with recoil techniques. 

In order to estimate the expected experimental sensitivity, we assume the same 
technique and the same performance in terms of efficiency and background rejection of 
the current BaBar search of Bd —>■ 77 [71] . The main source of background comes from 
hard photons from tt° or 77 decays, which are rejected requiring the invariant mass to be 
not within ±3<r of the nominal tt° or 77 mass. Other sources of background, as initial 
state radiation, are removed with cuts on topological variables. The efficiency on signal 
event is 10% and we expect 14 signal events and 20 background events in a sample of 
lab -1 of integrated data. We perform a set of toy Monte Carlo experiments, using tues, 
AE, and L2/L0 in a ML fit. The result in term of the integrated luminosity is shown 
in Fig. [TBI Considering that the measurement at lab -1 has already a 5a significance, 
we conclude that a super 5-Factory will be able to precisely measure the BR even at 
low statistic. Moreover, one will achieve a 7% statisticalprecision on the measurement 
at 30ab _1 , with a systematic error that will be smaller [] than 5%. The large amount 
of events available even with few ab _1 of integrated data will allow to measure also 
B s —>■ 77 CP asymmetry. This will be possible with good accuracy exploiting a tagging 
algorithm which determines whether the other B s meson in the event decayed as a B s or 
a B s (flavour tag) [72]. 

7 Current analyses of B° — > 7r°7r have a systematic error of 10%, generated by the presence of four 
photons in the final state and other sources of systematic, which have a statistical nature and can be 
strongly reduced at high integrated luminosity. 
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7 The impact on flavour physics: A possible future 
scenario 



In the previous sections we tried to give an idea of how rich is the set of physics mea- 
surements a B- Factory can perform while running at the T(5S). The processes described 
represent a minimum set of measurements demonstrating the possibility of accessing rates 
and weak phases of rare decays, i.e. the necessary ingredients to test the SM. 

Even with the limited amount of measurements we discussed, a facility like the one 
we are considering will have an impact on the knowledge of the flavour sector of NP. 
This might be clear enough from the previous chapters. In this section we make it more 
evident, showing how the unitarity triangle (UT) analysis will benefit from it. In order 
to do that, we consider two different scenarios: 

• a low- statistics case, assuming that after collecting 2 ab _1 of data at the T(4S"), (one 
of) the two existing 5-Factories will move to the T(5S'), collecting 1 ab _1 before 
LHCb will produce the first physics results. 

• a high- statistics case, in which a high luminosity 5-Factory will collect 75 ab _1 at 
the T(4S) and then 30 ab" 1 at the T(55) [13]. 

The low-statistics (high- statistics) scenario corresponds approximatively to fall 2009 (2015- 
2020). Taking into account the different time scales of the two projections, we consider 
two different sets of theoretical inputs: while for the low- statistics case we use the current 
values of the determinations of fs 3 ^Bb s , £, and B^, in the high- statistics case we assume 
that the increase of computation power will bring down the error of lattice QCD (LQCD) 
calculations to the percent level, while no improvement on the theory is considered. In- 
deed, since a constant progress in the calculation techniques has characterized the recent 
history of LQCD calculations, we can consider the estimates of the future errors relatively 
conservative. 

In Tab. EJ we quote the errors on the experimental and theoretical inputs currently 
used in the UT analysis, based on measurements performed at the T(4S"), together with 
6k from K—K mixing and Am s from B s -B s mixing. In Tab. H]we give the errors on the 
measurements performed at the T(5S), taken from the study presented in this paper. As 
mentioned in Sec. 15.2.11 we took into account an increase of the error of (3 S from J /ip<fi 
with tagged rates at positive and negative At, due to the presence of a two-fold ambiguity 
around /3 S ~ 0. 

7.1 The Unitarity Triangle in the Standard Model 

We show in Fig. [19] how the SM analysis of the UT will look like in the two scenarios 
described above. The two plots correspond to an overall error of o~p = 9% (cr^ = 4%) on 
p (fj) for the low-statistics scenario and to an overall error o~p = 1.4% (a^ = 0.8%) on p 
(fj) for the high- statistics one. For comparison, the current knowledge of the UT triangle 
gives o~p = 17% and = 5%. 
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Observable 


low- statistics 


high- statistics 






Error at 2 ab" 1 


Error at 75 ab -1 




v cb \ 


(exclusive) 


4% 


1.0% 




\v cb \ 


(inclusive) 


1% 


0.5% 




V u b: 


(exclusive) 


8% 


2.0% 




v ub 


(inclusive) 


8% 


2.0% 


BR(B -> tv) 


20% 


4% 




13% 


1% 






26% 


5% 






5% 10% 


1% 


sin(2/3) (J/ipK ) 


0.018 


0.005 






.Dif, combined) 


~ 6° 


1° 



a (combined) ~ 6° 1° 
2/3 + 7 IP^ttF, Z^ifgTrT) 20° 5° 



Table 3: Summary of the expected precision on the inputs of the UT analysis, in the low- 
statistics and high- statistics scenarios. The central values are taken in order to provide a 
perfect agreement among the constraints in the total fit. 



Observable 


low-statistics 


high- statistics 




Error at 1 ab^ 1 


Error at 30 ab" 1 


Ar 


0.16 ps" 1 


0.03 ps" 1 


r 


0.07 ps- 1 


0.01 ps" 1 


f3 s from angular analysis 


20° 


8° 




0.006 


0.004 


A C H 


0.004 


0.004 


P s from J/ip<j) At sign 


10° 


3° 


v td /v ts 


0.10 


0.031 


BR(B S -> fj+fj,-) ■ 10- s 


<10 at 90% prob. 


< 1.30 at 90% prob. 


BR(B S -> 77) 


38% 


7% 



Table 4: Expected errors for different observables considering an integrated luminosity of 
lab- 1 and SOab' 1 . 



7.2 The Unitarity Triangle beyond the Standard Model 

In a generic NP scenario the effect of AF = 2 NP contributions in B q -B q mixing can 
be parameterized in terms of two quantities, and <j>B , which relate the experimental 
observables to the SM quantities. They are defined as [73]: 



(B q \H™\ 


B q ) 


(B q \ 





—1 1 (91 1 
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p p 

Figure 19: Determination of p and fj in the Standard Model UT analysis for 2 ab _1 (left) 
and 75 atT 1 (right) projections of the errors as given in Tab. [3] 



where H^g is the effective Hamiltonian, is the Hamiltonian also including NP contri- 
butions, and AS M e 2i ^ M (A^e 2 ^™^^) is the SM (NP) amplitude. The two parameters 
C# and <pB q introduced in Eq. (I2T1) allow to relate the SM value of a certain observable 
to its measured value. For instance, the measured values of the size and phase of B^-B^ 
mixing are given by Amf xp = CgjAm^ and f3 EXP = j3 + <f>B d - Similar relations hold 
for the other observables of Tab. [3] Neglecting the case of NP contributions entering at 
tree-level processes, there are only two observables which do not depend on the presence 
of NP, namely \V ub /V cb \ and 7. E 

Since under these assumptions we can parameterize NP effects just with two real 
parameters, without assuming any specific model, it is possible to determine the allowed 
region for NP, fitting simultaneously CB d , 4>B d , p, and fj. Assuming the two scenarios of 
Tab. El we obtain the two plots given in Fig. [20] and the errors quoted in the first two 
raws of Tab. El The determination of Cs d and 4>B d in the two cases is given in Fig. [2TJ 

On top of this information, one can use the experimental information provided by the 
T(5S) run to constrain the values of Cb s and 4>b s - Cb s is directly related to the value of 
Am s , in analogy to the Bd case, while <pB a is defined in such a way that (3f xp = (3 S — 4>b s - 
The relation between the other constraints of Tab. H] and the NP parameters is more 
complicated. It has been shown that the semileptonic CP asymmetry [45] and the value 
of AT S /T S [47] are sensitive to NP contributions to AF = 2 effective Hamiltonian. We 
use the NLO expression of these observables in terms of the parameters Cb s and 4>b s [32]. 

The result of the combination of all these measurements is shown in the two plots 
of Fig. [22] and quantified in the bottom part of Tab. [5] On top of this analysis, one 
can further test the presence of NP using AF = 1 processes, as already discussed in 
Sec. 15.2.11 relating AF = 1 and AF = 2 NP effects by using a specific NP model (see for 
example [76]). 

8 In the case of 7, this is true up to a contribution in D-D mixing, which is estimated to be at the level 
of a few percent [74] and which can be taken into account following the approach suggested in ref. [75] . 
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p p 

Figure 20: Determination of p and fj in the generalized UT analysis, obtained using all 
the available constraints for 2 ab" 1 (left) and 75 atT 1 (right) projections of the errors as 
given in Tab. [3| 



Parameter 


low- statistics 


high- statistics 




Error 


Error 


P 


12% 


2.3% 


V 


8% 


1.8% 


C B d 


0.43 


0.032 


<pB d [°] 


1.7° 


0.44° 


Cb s 


0.33 


0.026 


<M°] 


7° 


1.9° 



Table 5: Determination of UT and NP parameters from the generalized UT analysis, in 
the two scenarios of Tab. For the B s sector, values from Tab. [7] are used. 

7.3 A comparison with the reach of the LHCb experiment 

Tab. [6] shows a comparison between the results presented in this paper and the expected 
sensitivity of LHCb for different benchmark measurement. The LHCb estimates are 
shown for an integrated luminosity of 2 fb^ 1 (one year of nominal data taking) and 
10 fb -1 (full data taking), using the errors on lattice predictions expected on 2009 and 
2015 respectively. The two scenarios described in the previous sections are used for 
a Z?-Factory running at the T(5S) resonance. The measurement of Am s at LHCb is 
supposed to be dominated by a systematic error of 0.09 ps~ Y [39]. The relative error on 
Bfj — ► K* 0/ -f is estimated to be ~ 5% (dominated by the systematics on the branching 
ratio's normalization). The relative error on B s — > K* 0/ y is estimated to be ~ 20% in 
the 2 scenario and ~ 10% in the 10 fb^ 1 scenario [77,78]. In the A S SL and Ach 
measurements both statistical and systematic uncertainties are shown. 

In doing this comparison, the reader should remember that a reliable estimate of im- 
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Figure 21: Determination of Cb a and 4>B d in the generalized UT analysis for 2 ab -1 (left) 
and 75 ab -1 (right) projections of the errors given in Tab.\^ combined to the low-statistics 
(left) and high- statistics (right) scenarios of Tab. 




Figure 22: Determination of Cb 3 and <pB 3 in the generalized UT analysis for 1 ab -1 (left) 
and 30 ab -1 (right) projections of the errors for the input of Tab. [^J In the left (right) 
plot, 2 ab -1 (75 ab^ 1 ) inputs of Tab. [3] are used for the other determinations. 



portant systematic effects is not available for several LHCb measurements, while these 
uncertainties are typically well known and under control at the I?-Factories. Moreover, we 
want to stress here that we presented only a minimal subsample of the possible measure- 
ments that can be performed at a S-Factory running at the T(5S), aiming to demonstrate 
the experimental feasibility of the corresponding analysis techniques. In the channels we 
chose, our estimates are mostly not competitive with respect to LHCb, but through them 
we proved the possibility to widely apply these techniques, in order to investigate several 
channels that cannot be studied in the LHC environment. 
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LHCb 


T(5S) 


Observable 




10 /& 1 


lafr -1 




|Kd/Ks| from Am s 


n m n 


002 






AT/T 


fl DDQ9 


004 


1 9 


09 


3 S from angular analysis 


0.66° 


0.29° 


20° 


8° 








±0.006 ±0.004 


±0.001 ±0.004 


A C H 






±0.0015 ±0.004 


±0.0003 ±0.004 


8 S from J/ip(f) At sign 






20° 


8° 


BR(B S -> yu/i) 


1.2 • 1(T 9 


0.7- 10- 9 


< 10- 7 


< 1.30 • 10- 8 


^/14 s from radiative decays 


0.03 


0.015 


0.10 


0.031 


BR(B S -> 77) 






38% 


7% 



Table 6: Expected errors for different observables at LHCb [39, 11, 18] and at a B-Factory 
running at the T(5S) resonance. 



8 Conclusions 

In this paper we discussed the physics potential of a 5-Factory running at the T(5S). 
Assuming the performance of the detectors currently taking data at the two 5-Factories, 
we proved that it is possible to obtain a pure sample of B s mesons and at the same time 
to cover the physics program of a traditional 5-Factory. We discussed the detector limi- 
tations in studying B s -B s oscillations and we showed that, even if time-dependent mea- 
surements of B s are not accessible with the current vertex resolution, the same physics 
information can be obtained using time-integrated measurements. We also illustrated 
some benchmark analyses of rare B s decays, interesting to test the SM and obtain infor- 
mation on the flavour structure of the NP. With this limited set of examples, we proved 
that this facility can perform similar measurements using alternative techniques with re- 
spect to hadron collider experiments. On the other hand, thanks to the clean environment 
and the progress made by BaBar and Belle in the recent past, this facility allows to ex- 
tend the LHCb physics program to a larger set of interesting channels which, not having 
prompt charged tracks and/or being characterized by the presence of photons in the final 
state, are difficult to study at a hadron collider. Examples of this kind are B s — > 77 
and B s — > K°K°. To conclude, we extrapolated the impact of the B s physics program 
on the unitarity triangle analysis beyond the standard model, considering the possibility 
of lab -1 collected by the existing 5-Factories after accumulating 2ab -1 at the T(4S), or 
30ab -1 collected by a high luminosity _B-Factory after collecting 75ab _1 at the T(AS). We 
stress that the experimental precision on these observables is not necessarily better than 
what LHCb can achieve [79]. Nevertheless, this limited set of examples demonstrates 
that, even using different experimental techniques, the proposed machine can measure 
absolute values and (strong and weak) phases of decay amplitudes, even without directly 
accessing the time-dependent structure of B s -B s oscillation. Furthermore, thanks to the 
very clean environment, it is possible to apply the discussed experimental techniques to all 
B decays, without any restriction related to the need of prompt tracks in the event. We 
do not cover all the examples of such additional potentiality with respect to hadron col- 
lider experiments, but only outline the cases of the measurements of the rate of B s — > 77 
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decays and of the weak phase in B s — > K°K° decays. 

Concluding, we would like to stress that, unlike the case of B d and charged B decays, 
the B s sector is presently largely unstudied. Based on what we discussed in this paper, we 
think that it is particularly interesting to consider the possibility of studying B s decays 
with BaBar and Belle before LHCb starts. This would maximize the physics impact of 
the luminosity that can be collected (less than 2ab _1 ), while a super B- Factory would be 
needed to obtain the required luminosity during and after LHC data taking. 

9 Acknowledgments 

We would like to thank G. Isidori and G. Martinelli for their help in writing Sec. I4.1[ 
V. Lubicz for providing us the extrapolation of lattice QCD errors, R. Faccini for useful 
discussions and suggestions and for the feed-back provided during all the development of 
this work and E. Emili and D. Granata for their contribution in our early studies for this 
paper. We would also like to thank T. Browder and M. Hazumi for the interest shown 
for this work and for inviting us to take part at the BNMS2006 workshop series, which 
allowed us to compare our studies to the activity of the Belle community. This work 
has been supported in part by the France- Japon Particle Physics Laboratory FJ-PPL, 
the EU networks "The quest for unification" under the contract MRTN-CT-2004-503369, 
"FLAVIAnet" under the contract MRTN-CT-2006-035482, and "HEPTOOLS" under the 
contract MRTN-CT-2006-035505. 

References 

[1] N. Cabibbo, Phys. Rev. Lett. 10 (1963) 531. M. Kobayashi and T. Maskawa, Prog. 
Theor. Phys. 49 (1973) 652. 

[2] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 87 (2001) 091801 
[arXiv:hep-ex/0107013]. K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 87 
(2001) 091802 [arXiv:hep-ex/0107061]. B. Aubert et al. [BABAR Collaboration], 
arXiv:hep-ex/0607107. K. F. Chen et al. [Belle Collaboration], arXiv:hep-ex/0608039. 

[3] M. Ciuchini, E. Franco, G. Martinelli, L. Reina and L. Silvestrini, Z. Phys. C 68 
(1995) 239 [arXiv:hep-ph/9501265]. 

[4] P. Paganini, F. Parodi, P. Roudeau and A. Stocchi, Phys. Scripta 58 (1998) 556 
[arXiv:hep-ph/9711261]. 

[5] M. Bona et al. [UTfit Collaboration], JHEP 0610 (2006) 081 [arXiv:hep-ph/0606167]. 
Updated results are available on http://www.utfit.org. For a different approach, 
producing similar results, see J. Charles et al. [CKMfitter Group], Eur. Phys. J. C 
41 (2005) 1 [arXiv:hep-ph/0406184]. 

[6] B. Aubert et al. [BABAR Collaboration], Nucl. Instrum. Meth. A 479 (2002) 1 
[arXiv:hep-ex/0105044] . 

[7] A. Bondar et al. [BELLE Collaboration], Nucl. Instrum. Meth. A 408 (1998) 64. 



43 



[8] F. Abe et al. [CDF Collaboration], Nucl. Instrum. Meth. A 271, (1988) 387. J. Cran- 
shaw [CDF Collaboration], Nucl. Instrum. Meth. A 462, (2001) 170. 

[9] S. Abachi et al. [D0 Collaboration], Nucl. Instrum. Meth. A 338, (1994) 185. 
S. Abachi et al. [D0 Collaboration], Nucl. Instrum. Meth. A 408, (1998) 103. 

[10] V. M. Abazov et al. [D0 Collaboration], Phys. Rev. Lett. 97 (2006) 021802 
[arXiv:hep-ex/0603029]. A. Abulencia [CDF - Run II Collaboration], Phys. Rev. Lett. 
97 (2006) 062003 [arXiv:hep-ex/0606027]. 

[11] O. Schneider et al. [LHCb Collaboration], Nucl. Instrum. Meth. A 446, 213 (2000). 

[12] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 71 (2005) 032005 [arXiv:hep- 
ex/0411016]. R. Itoh et al. [Belle Collaboration], Phys. Rev. Lett. 95 (2005) 091601 
[arXiv:hep-ex/0504030]. 

[13] M. Bona, "SuperB: A high-luminosity asymmetric e+ e- super flavor factory. Con- 
ceptual design report", URL http://www.pi . infn. it/SuperB/?q=CDR. 

[14] W. M. W. Yao et al. [Particle Data Group], J. Phys. G 33 (2006) 1. 

[15] G. Huang et al. [CLEO Collaboration], arXiv:hep-ex/0607080. 

[16] D. Besson et al. [CLEO Collaboration], Phys. Rev. Lett. 54 (1985) 381. D. M. J. Love- 
lock et al, Phys. Rev. Lett. 54 (1985) 377. 

[17] A. Drutskoy, for Belle Collaboration arXiv:hep-ex/0605110. 

[18] See for example N. A. Tornqvist, Phys. Rev. Lett. 53 (1984) 878. 

[19] J. Ocariz et al, Babar Analysis Document #346 

[20] R. Fisher, Ann. Eugenics, vol. 7, pp. 179-188, 1936. 

[21] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 71 (2005) 111102 [arXiv:hep- 
ex/0503011]. 

[22] See for example B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 71 (2005) 
091102 [arXiv:hep-ex/0502019]. 

[23] S. J. Brodsky and J. R. Hiller, Phys. Rev. D 46 (1992) 2141. 

[24] F. Cardarelli, I. L. Grach, I. M. Narodetsky, G. Salme and S. Simula, Phys. Lett. 
B 349 (1995) 393 [arXiv:hep-ph/9502360]. D. Melikhov and S. Simula, Phys. Rev. 
D 65 (2002) 094043 [arXiv:hep-ph/0112044]. B. L. G. Bakker, H. M. G. Choi and 
C. R. G. Ji, Phys. Rev. D 67 (2003) 113007 [arXiv:hep-ph/0303002]. 

[25] G. C. Branco, L. Lavoura, and J. P. Silva, "CP Violation", Oxford University Press 
(1999). 

[26] M. Gronau and D. London, Phys. Rev. Lett. 65 (1990) 3381. 



44 



[27] M. Bona et al. [UTfit Collaboration], arXiv:hep-ph/0701204. 

[28] B. Aubert et al. [BABAR Collaboration], arXiv:hep-ex/0607106. 

[29] P. Oddone, Annals N. Y. Acad. Sci. 578 (1989) 237. 

[30] B. Aubert et al. [BABAR Collaboration], arXiv:hep-ex/0203035. 

[31] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 66 (2002) 032003 [arXiv:hep- 
ex/0201020]. 

[32] M. Bona et al. [UTfit Collaboration], JHEP 0603 (2006) 080 [arXiv:hep-ph/0509219]. 
M. Bona et al. [UTfit Collaboration], Phys. Rev. Lett. 97 (2006) 151803 [arXiv:hep- 
ph/0605213]. 

[33] D. Acosta et al. [CDF Collaboration], Phys. Rev. Lett. 94 (2005) 101803 [arXiv:hep- 
ex/0412057]. 

[34] H. Albrecht et al. [ARGUS Collaboration], Z. Phys. C 48 (1990) 543. 

[35] G. Batignani et al, contribution to the 15th International Work- 
shop on Vertex Detectors (Vertex06), available at the URL 
http : // indico . cern . ch/materialDisplay . py?contribId=31&amp ; sessionld=10 
&amp ; mater ialld=paper&amp ; conf Id=1151. 

[36] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 93 (2004) 131805 
[arXiv:hep-ex/0403001]. 

[37] V. M. Abazov et al. [D0 Collaboration], arXiv:hep-ex/0701012. 

[38] S. Amato et al. [LHCb Collaboration], CERN-LHCC-98-4 

[39] M. H. Schune, talk given at the 5 th Super B Workshop Paris, May 9-11, 2007. 

[40] M. Ciuchini, M. Pierini and L. Silvestrini, Phys. Rev. Lett. 95 (2005) 221804 
[arXiv:hep-ph/0507290]. 

[41] M. Ciuchini, M. Pierini and L. Silvestrini, arXiv:hep-ph/0703137. 

[42] A. J. Buras and L. Silvestrini, Nucl. Phys. B 569 (2000) 3 [arXiv:hep-ph/9812392]. 

[43] N. Neri and M. Pierini, study for the superB factory, presented at the first super B 
workshop, Frascati, Italy 11-12 November 2005 

[44] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 97 (2006) 171805 
[arXiv:hep-ex/0608036]. 

[45] S. Laplace, Z. Ligeti, Y. Nir and G. Perez, Phys. Rev. D 65 (2002) 094040 [arXiv:hep- 
ph/0202010]. 

[46] D. Acosta et al. [CDF Collaboration], "Analysis of decay-time dependence of angular 
distributions in B® — > J /ip<fi and B® — > J/if)K*° decays and measurement of the Phys. 
Rev. Lett. 94 (2005) 101803 [arXiv:hep-ex/0412057]. 



45 



[47] A. S. Dighe, I. Dunietz and R. Fleischer, Eur. Phys. J. C 6 (1999) 647 [arXiv:hep- 
ph/9804253]. I. Dunietz, R. Fleischer and U. Nierste, Phys. Rev. D 63 (2001) 114015 
[arXiv:hep-ph/0012219]. 

[48] B. Aubert et al. [BABAR Collaboration], arXiv:hep-ex/0607067. 

[49] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 96 (2006) 251802 
[arXiv:hep-ex/0603053]. 

[50] B. Grinstein and R. F. Lebed, Phys. Rev. D 53 (1996) 6344 [arXiv:hep-ph/9602218]. 

[51] Y. Grossman, Z. Ligeti, Y. Nir and H. Quinn, Phys. Rev. D 68 (2003) 015004 
[arXiv:hep-ph/0303171]. 

[52] M. Ciuchini, M. Pierini and L. Silvestrini, Phys. Lett. B 645 (2007) 201 [arXiv:hep- 
ph/0602207]. 

[53] A. Pomarol and D. Tommasini, Nucl. Phys. B 466 (1996) 3 [arXiv:hep-ph/9507462]. 
R. Barbieri, G. R. Dvali and L. J. Hall, Phys. Lett. B 377 (1996) 76 [arXiv:hep- 
ph/9512388]. R. Barbieri, L. J. Hall and A. Romanino, Phys. Lett. B 401 (1997) 47 
[arXiv:hep-ph/9702315]. 

[54] M. Ciuchini et al, JHEP 0107 (2001) 013 [arXiv:hep-ph/0012308]. M. Bona et al. 
[UTfit Collaboration], JHEP 0507 (2005) 028 [arXiv:hep-ph/0501199]. 

[55] M. Beneke, T. Feldmann and D. Seidel, Nucl. Phys. B 612 (2001) 25 [arXiv:hep- 
ph/0106067]. S. W. Bosch and G. Buchalla, Nucl. Phys. B 621 (2002) 459 [arXiv:hep- 
ph/0106081]. A. Ali, E. Lunghi and A. Y. Parkhomenko, Phys. Lett. B 595, (2004) 
323 [arXiv:hep-ph/0405075]. S. W. Bosch and G. Buchalla, JHEP 0501 (2005) 035 
[arXiv:hep-ph/0408231]. M. Beneke, T. Feldmann and D. Seidel, Eur. Phys. J. C 41 
(2005) 173 [arXiv:hep-ph/0412400]. 

[56] M. Nakao et al. [BELLE Collaboration], Phys. Rev. D 69 (2004) 112001 [arXiv:hep- 
ex/0402042]. B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 70 (2004) 
112006 [arXiv:hep-ex/0407003]. 

[57] T. E. Coan et al. [CLEO Collaboration], Phys. Rev. Lett. 84 (2000) 5283 
[arXiv:hep-ex/99 12057]. K. Abe et al, Phys. Rev. Lett. 96 (2006) 221601 [arXiv:hep- 
ex/0506079]. B. Aubert [BABAR Collaboration], Phys. Rev. D 66 (2002) 010001 
arXiv:hep-ex/0612017. 

[58] P. Ball and R. Zwicky, JHEP 0604 (2006) 046 [arXiv:hep-ph/0603232]. 

[59] A. J. Buras, M. Spranger and A. Weiler, Nucl. Phys. B 660 (2003) 225 [arXiv:hep- 
ph/0212143]. 

[60] T. Inami and C.S. Lim, Prog. Theor. Phys. 65 (1981) 297; ibid. 65 (1981) 1772. 
[61] M. Blanke, A. J. Buras, D. Guadagnoli and C. Tarantino, arXiv:hep-ph/0604057. 



46 



[62] C. Bobeth, M. Bona, A. J. Buras, T. Ewerth, M. Pierini, L. Silvestrini and A. Weiler, 
Nucl. Phys. B 726 (2005) 252 [arXiv:hep-ph/0505110]. 

[63] G. Isidori and A. Retico, JHEP 0111 (2001) 001 [arXiv:hep-ph/0110121]. 

[64] T. Bergfeld et al. [CLEO Collaboration], Phys. Rev. D 62 (2000) 091102 [arXiv:hep- 
ex/0007042]. M. C. Chang et al. [BELLE Collaboration], Phys. Rev. D 68 (2003) 
111101 [arXiv:hep-ex/0309069]. B. Aubert et al. [BABAR Collaboration], Phys. Rev. 
Lett. 94 (2005) 221803 [arXiv:hep-ex/0408096]. 

[65] CDF Collaboration, note #8176. 

[66] M. Ciuchini, E. Franco, A. Masiero and L. Silvestrini, Phys. Rev. D 67 (2003) 075016 
[Erratum-ibid. D 68 (2003) 079901] [arXiv:hep-ph/0212397]. 

[67] L. Reina, G. Ricciardi and A. Soni, Phys. Rev. D 56 (1997) 5805 [arXiv:hep- 
ph/9706253]. 

[68] A. Gemintern, S. Bar-Shalom and G. Eilam, Phys. Rev. D 70 (2004) 035008 
[arXiv:hep-ph/0404152]. 

[69] S. Bertolini and J. Matias, Phys. Rev. D 57 (1998) 4197 [arXiv:hep-ph/9709330]. 

[70] G. L. Lin, J. Liu and Y. P. Yao, Phys. Rev. Lett. 64, 1498 (1990). G. L. Lin, J. Liu 
and Y. P. Yao, Phys. Rev. D 42 (1990) 2314. 

[71] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 87 (2001) 241803 
[arXiv:hep-ex/0107068]. 

[72] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 92 (2004) 201802 
[arXiv:hep-ex/0312055]. 

[73] J. M. Soares and L. Wolfenstein, Phys. Rev. D 47 (1993) 1021. N. G. Deshpande, 
B. Dutta and S. Oh, Phys. Rev. Lett. 77 (1996) 4499 [arXiv:hep-ph/9608231]. 
J. P. Silva and L. Wolfenstein, Phys. Rev. D 55 (1997) 5331 [arXiv:hep-ph/9610208]. 
A. G. Cohen, D. B. Kaplan, F. Lepeintre and A. E. Nelson, Phys. Rev. Lett. 78 
(1997) 2300 [arXiv:hep-ph/9610252]. Y. Grossman, Y. Nir and M. P. Worah, Phys. 
Lett. B 407 (1997) 307 [arXiv:hep-ph/9704287]. 

[74] A. Giri, Y. Grossman, A. Soffer and J. Zupan, Phys. Rev. D 68 (2003) 054018 
[arXiv:hep-ph/0303187]. 

[75] A. Amorim, M. G. Santos and J. P. Silva, Phys. Rev. D 59 (1999) 056001 [arXiv:hep- 
ph/9807364]. 

[76] L. Silvestrini, arXiv:0705.1624 [hep-ph]. 

[77] J. Dickens, talk given at the CKM 2006 Workshop, Nagoya, December 12-16, 2006. 
[78] V. Vagnoni, private communication. 



47 



[79] O. Schneider, talk given at the LHCb Upgrade Workshop Ed- 
inburgh, January 11-12, 2006. Slides available at the URL 
http : //lphe . epf 1 . ch/publications/2007/0Schneider_Edinburgh . ppt 



48 



